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Abstract 



Particle Image Velocimetry (PIV) is a non-intrusive technique for simultaneously 
measuring the velocities at many points in a fluid flow. The fluid is seeded with 
tracer particles and the region under investigation is illuminated. An image of the 
illuminated region is captured and then, a short time period later, a second image is 
taken. Suitable analysis of these images yields an instantaneous velocity vector map. 

Until recently, restrictions in the rate at which images could be captured have 
limited the PPV technique to the analysis of slow flows. However, advances in camera 
technology have now opened up the possibility of using PPV in the analysis of faster 
flows. Indeed, image capture rates are now fast enough to enable two images to be 
captured during a fraction of an acoustic cycle, indicating the potential for using PIV 
to analyse sound fields. 

In this thesis, after some aspects of sound field theory have been outlined and 
following a discussion of the theory of PIV, the development of experimental PIV 
apparatus for measuring sound fields is described. Measurements of the temporal 
variation in the velocities of particles within some common sound fields are presented. 
In particular, the passage of an acoustic pulse is monitored and the sinusoidal motion 
of particles in a resonating tube is recorded yielding the corresponding standing wave 
pattern. Finally, the main limitations of the PIV technique when applied to acoustic 
fields are discussed. 
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Introduction 



After the first acoustic velocity measurements were carried out by 
Rayleigh at the end of the 19 th century, there followed a long period 
during which there were only small advances in the measurement of ve- 
locities in sound fields. Indeed, it was nearly a century later, with the employ- 
ment of optical methods, that acoustic velocimetry took its next leap forward. 

Particle Image Velocimetry (PIV) has become a widely used optical method 
for measuring fluid flow velocities. It has the unique features of being non- 
intrusive and whole-field; that is, simultaneously measuring the velocities at 
many points in a fluid flow without any disturbance to it. The technique in- 
volves seeding the fluid with tracer particles and illuminating the region under 
investigation. Two or more images of the illuminated region are then captured 
in rapid succession. Finally, the images are analysed to extract the velocity values. 

Until recently, restrictions in the rate at which images could be captured have 
limited the PIV technique to the analysis of slow flows. However, advances in 
camera technology have now opened up the possibility of using PIV in the analy- 
sis of faster flows, in particular for the measurement of sound fields. 



1.1 The Aims 

The main goals of the present work are: 

1. To investigate the possibility of employing PIV to study high-speed acoustic 
motions, exploring the limits under which the principles underpinning PIV 
remain valid for acoustic measurements. 
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2. To develop a PIV set-up which is capable of measuring acoustic particle ve- 
locities. That is, a measurement system that is able to accurately record the 
oscillatory motion of particles in sound fields. 

3. To employ the extended PIV set-up in the spatial and temporal measurement 
of acoustic velocity. 

1.2 Outline of Thesis 

Chapter 2 presents a brief theory of sound fields. In particular, the equations of 
resonant tubes, which are later studied experimentally, are introduced. 

Chapter 3 is an in-depth review of the theory and development of Particle 
Image -Velocimetry The chapter covers the structure of PIV systems, describing 
in detail each component of the technique. It then presents the theory behind 
the most common method of PIV data analysis. Throughout, specific attention 
is paid to the digital form of the technique. 

Based on the previous chapters, Chapter 4 discusses the potential problems 
associated with using PIV for sound field measurements. The assumptions upon 
which PIV are based are examined to verify that they remain valid when measur- 
ing sound fields. In particular, the level of intrusion caused by the measurements 
and the degree to which the results truly reflect the actual acoustic velocities are 
investigated. 

Chapter 5 describes the experimental configuration used in this study pre- 
senting the design of the data acquisition and data analysis systems. 

Chapter 6 examines some applications of acoustic PIV. It presents the results 
of the measurement of some common sound fields. In particular, the passage of 
an acoustic pulse is monitored and the sinusoidal motion of particles in resonat- 
ing tube is recorded yielding the corresponding standing wave pattern. Finally, 
it discusses some unexpected results encountered during the acoustic PIV mea- 
surement of oscillatory motions. 

The thesis concludes with Chapter 7 which contains a summary of the main 
findings of the research and some ideas for future work. 
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Sound Fields 



All science is either 

physics or 

stamp collecting. 



Ernest Rutherford 



T 



his chapter covers those parts of sound field theory which are related to 
this study. A detailed development of the theory presented in this chap- 
ter can be found in most standard acoustics and fluid dynamics text- 



books ( |Meyer and Neumann||i972[ |Morse and Ingard||i968| |Kinsler et al, 
I2000 



Skudrzyk 1971; Landau and Lifshitz 1987 1 and consequently only the 



most pertinent results are outlined here. 



2.1 Terminology 



2.1.1 



Fluid Flows 



Throughout this study, it is considered that the fluid under investigation con- 
sists of fluid elements, which are volumes small enough that all the sound field 
parameters can be assumed to be constant within them. The position of a fluid 
element, in the Cartesian co-ordinate system, is denoted by 



[x,y,z] = x$ + -y$ + zz, 



(z.1.1) 



where x, y and z are the unit vectors. When at rest, i.e. when there is no force 
applied to any of its elements, the fluid has an equilibrium pressure of To and 
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equilibrium density of po. 

When the effect of energy dissipation is considered, the fluid is said to be 
real or viscous. Energy dissipation in fluids is an irreversible process which is 
caused by viscosity, which is the internal friction between molecules of the fluid, 
or thermal conduction or both. In the idealised case, where all of these effects are 
negligible, the fluid is referred to as ideal or inviscid. 

Fluids are also classified according to their compressibility. In a compressible 
fluid the equilibrium density is a function of space and time while it is a constant 
value in an incompressible fluid. 

2.1.2 Acoustical Parameters 

In the presence of forces generated by a sound field, the instantaneous properties 
of the fluid are said to be acoustical in nature. Under this condition, a fluid 
element has an instantaneous displacement from its equilibrium position denoted 
by 

The velocity associated with this instantaneous displacement is known as the 
acoustic velocity and is denoted 

v"= [v x ,v y ,v z ) = 3 t £,. (2.1.2) 

The acoustic pressure p is the difference between the the instantaneous pressure 7 
and the equilibrium pressure: 

V = 7-? . (2.1.3) 

The instantaneous density p defines the condensation of the fluid as 

P - Po 1 X 

s = . (2.1.4) 

Po 

2.2 Fundamental Equations 

A sound field in a fluid can be fully described by the corresponding acoustic 
velocity V and two thermodynamic quantities, which are usually acoustic pres- 



2.2. Fundamental Equations 



sure CP and condensation s. These parameters are related by three fundamental 
equations presented in this section. 

2.2.1 Equation of Motion 

The equivalent of Newton's equation of motion in an incompressible viscous 
fluid is the Navier-Stokes equation which in its general form may be written as 

p (a t v +(v.V) : v) =-Vp+riV 2 ^+ (-ti + tibJ V(V.v). (2.2.1) 

where n is the shear or dynamic viscosity coefficient of the fluid and t)b is the bulk 
viscosity coefficient of the fluid. 

For an incompressible fluid, when V. v" = and the viscosity is zero, the 
Navier-Stokes equation is reduced to Euler's equation of inviscid motion given 

by 

p ( d t v"+ ( v\V) vj = — Vp. (2.2.2) 

For steady flow, Euler's equation simplifies to Fourier's law, which is expressed 
as 

p( v\V) v" = — Vp. (2.-2.3) 

In the linear acoustics approximation, the second order term in Euler's equa- 
tion is negligible and the equation is reduced to 

p3 t v" = -Vp, (2.2.4) 

which is known as the linear Euler equation. 

2.2.2 Equation of Continuity 

Based on the law of conservation of mass, the equation of continuity of a fluid can 
be derived as 

9 tP + V.(pv) = 0, (2.2.5) 

which relates the motion of the fluid to its compression or expansion. In the 
linear acoustic estimation, the linear equation of continuity is given by 

a t s + V.\T = 0. (2.2.6) 
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2.2.3 Equation of State 



When the host fluid is a perfect gas, in most cases, the change of all parameters 
of the sound field is much faster than the change of temperature. As a result, 
the temperature can be assumed to be a constant implying that the host fluid is 
adiabatic. This requires the condition 

where y is the ratio of specific heats given by 

y = —, (2.2.8) 

c v 

with cp and cy representing the specific heat at constant pressure and specific heat 
at constant volume, respectively. 

The equation of state for an adabatic gas is given by 

p = c 2 (p-p ), (2.2.9) 

where c is the thermodynamic speed of sound given by 

c = — . (2.2.10) 

V Po / 



In the linear approximation, when the condensation is much smaller than 
unity, the linear equation of state can be found from Taylor's expansion of the 
instantaneous pressure around its equilibrium value. This yields 

p = 2 s, (2.2.L1) 

where S is the adiabatic bulk mudulus given by 

S = Po(3 p ?) Po . (2.2.12) 

8 



2j. Linear Sound Waves 

2.3 Linear Sound Waves 

2.3.1 Propagating Waves 

Combining the linear forms of the fundamental equations for an inviscid fluid 



yields the linear lossless wave equation of propagation of sound as (Kinsler et al. 
\1000) 

D 2 p(x,t)=0, (2.3.1) 

where D 2 denotes the d'Alembertian operator defined as 

□ 2 = v 2 -c- 2 a t 2 , (2.3.2) 

with c being the speed of propagation of sound in the fluid and V 2 denoting the 
Laplacian operator. 

The general solution of Equation (2.3.1]) is a harmonic plane wave. That is, a 



wave where, at any point on a plane perpendicular to the direction of propaga- 
tion, the properties of the sound field are the same. 

In potential or irrotational fluids, i.e. when V x v = 0, the acoustic veloc- 
ity may be presented in terms of the gradient of a scalar, known as the velocity 
potential, by 1 

v"(x",t) = V*(x,t). (2.3.3) 



Using the linear Euler equation (2. 2. 4)), the pressure can be expressed in terms 



of the temporal differentiation of the velocity potential: 

p(x\t) =-p d t <k[x,t). (2.34) 

Practically, when calculating properties of sound fields, the aim is to obtain the 
velocity potential which in turn yields the other parameters. It can be shown 
that, with the same linear approximation, the velocity potential satisfies the wave 
equation as 

□ 2 c|)(x,t)=0, (2.3.5) 

In this study, the plane waves of interest are one-dimensional, satisfying the 



'Sometimes the velocity is expressed, conventionally, by the gradient of the velocity potential 
with a minus sign, i.e. V [x, t) = — Vc|>(x, t). 
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reduced linear wave equation 

(d 2 x -c- 2 d t 2 )<b(x,t)=0, (2.3.6) 

which has real harmonic solutions in terms of the trigonometric functions. It is 
convenient to express these solutions in the complex exponential form and then 



use the real part. 2 With this in mind, the solution of Equation (2.3.6! can be 
written as 

*(x,t)=JK[<Mx)e- UDt ], (^3-7) 

where to is the angular frequency of the wave and the velocity potential cj) x is a 
function only of space fulfilling the Helmholtz equation 

(V 2 + k 2 )cM5?)=0, (2.3-8) 

with k being the wave number related to the wavelength A by 

k= — = — . (2.3.9) 

c A 



Solving the Helmholtz equation (2.3.8 1, the final form of the velocity potential 
may be given by 

4>(x, t) = m [<t> e x{kx - wt] ] , (2.3.10) 

where cj)o is the amplitude. Equation ( 2.3.L0I describes a monochromatic (i.e. with 



constant frequency cu) travelling wave propagating in an infinite medium in the 
x direction and with a speed of tu/k. 

The linear lossy wave equation for a viscous fluid may be obtained by combin- 
ing the linear equations of continuity and state with the Navier- Stokes equation 
for incompressible fluid to give 

D 2 c|)(x ) t)=0, (2.3.11) 

where □£ stands for the lossy d'Alembertian operator defined as 

□ 2 = U + T f a t )v 2 -c- 2 d t 2 , (2.3.12) 



2 Usually, when only linear operations are involved, the real part sign is omitted during cal- 
culations but appears in the final results. 
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with Tf denoting the relaxation time of the viscous fluid element. 

It can be shown (Kinsler et al.||20oo §8.2) that the general solution of the 



lossy wave equation ( 2.3.11I is the same as the lossless version but with a complex 
wave number 

k = k + x<x, (2..3.13) 

where the boldface denotes a complex variable and a is the absorption coefficient. 
Thus, one-dimensional wave propagation in a real fluid can be described by 

<P(x,t) = <$>oe x{kx - a,t] . (2.3.14) 



2.3.2 Standing Waves 

Many observed acoustical phenomena take place in finite media where the sound 
waves meet the boundaries of the finite media. In such a case, based on the prop- 
erties of the boundaries, the incident waves convert to reflected and/or transmit- 
ted waves. 

Sound fields in finite one-dimensional media can be described by acoustic 



transmission line theory. With reference to Equation (2.3.10 1, a lossy transmission 
line consists of an incident wave 

^{xA) = ^e x(+ ^- wX \ (2.3.15) 

and its reflection from a real boundary 

^{x,t) = ^t^ [ ^~ ( " t \ (2.3.16) 

where it is assumed that the boundary wall is rigid or acoustically hard, that is, the 
incident wave is almost completely reflected. The linear superposition of these 
waves at each point inside the finite medium and at any instant yields 

<j>(x,t) = (|j- > (x,t) + <|) < -(x,t). (2.3-17) 



Applying the boundary conditions for a finite medium, Equation (2.3.17) can 
only be satisfied for a limited number of frequencies known as the natural or 
characteristic frequencies of that finite medium. Moreover, the waves described 
by the solutions of this equation do not propagate in space and are said to be 
stationary or standing waves. 
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2.3.3 Impedance 

The acoustic impedance acting on a surface of area S is defined as the quotient of 
the complex acoustic pressure and the complex volume velocity associated with 
that area: 

Z[x,t) = P^%. (2.3.18) 

Sv(x,t) 

The specific acoustic impedance is, in general, given by the quotient of the acoustic 
pressure and velocity as 

Z s (x,t) = *M (2.3.19) 

v(x,t) 

In the special case, where the acoustic pressure and velocity correspond to a plane 
wave, the specific acoustic impedance is 

Z = Poc. (2.3.20) 

which is known as the characteristic acoustic impedance or the wave impedance of 
the medium. 



2.4 Theory of the Resonance Tube with Rigid Ends 

Consider a tube 3 of length L which is driven at one end (x = 0) by a sound 
generating device, eg a piston or a loudspeaker, and closed at the other end (x = 
L). Both ends of the tube are assumed to be rigidly terminated. Let the sound 
generating device produce a monochromatic wave with frequency a> shown by 

x'(t) = Isin(tut), (2.4.1) 

where I is the displacement amplitude of the sound generating device. 

In the steady state, the velocity potential created by the sound source inside 
the tube, as a finite medium, can be described by the standing wave Equation 



( I2.3.17P as 

* (x, t) = c|V e t(+1cx -' ut) + 4C e x{ - kx - wt] . (2.4.2) 



3 Other commonly used synonyms of 'tube' in the literature are duct, pipe or wave-guide. 



12 



2.4- Theory of the Resonance Tube with Rigid Ends 

The 'exact' boundary conditions of this system may be given by 

v(x',t) = 3 x (|)(x,t) \ x=x > = 3 t x'(t), (2.4-3) 

for the sound generating side and 

v(L,t) = d x (\)(x,t) U=l = 0, (2.44) 

for the rigid termination side of the tube. 



2.4.1 Assumptions 



Before solving the standing wave Equation (2.4.2 1 some assumptions must be 
taken into account. 

First, it is supposed that the wave equation is linearised. However, the bound- 



ary condition Equation (2.4.3 ' leads the constant velocity amplitudes v^ and v^ 
to be a function of time and, consequently, generating non-linear solutions. This 
may be overcome by assuming that the displacement amplitude of the sound 



generating device is negligible comparing to the length of the tube (Temkin 198 
§3.6) or 

T < 1. ( 2 -4-5) 



which causes the boundary condition Equation (2.4.3 1 to De replaced by 



v(0,t) = 9 x cMx,t) | x=0 = 0. 



(2.4.6) 



Hence, Equations (2.4. 5 1 and (2.4.6 1 are required to be met for a linear regime. 

The second assumption is important in order to have a plane and longitudi- 
nal sound wave inside the tube. Of these criteria, the first one can be fulfilled 
when the length of the tube is much larger than its transverse dimensions or, in 
other words, if 

L > VS. (2.4.7) 

where S denotes the transverse cross section of the tube. Apparently, the wave 
created by the sound generation device can be purely longitudinal when the de- 
vice motion is perpendicular to the axis of the tube. However, in practice this 
may not happen perfectly. It can be shown that there exists a cut-off frequency 
proportional to l/VS under which, even when the device motion normal is not 
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absolutely axial, any transverse waves are dumped quickly (Temkin 1981, §3.5). 
This criterion may be written as 



to < 



I 

7s" 



(2.4.8) 



Finally it is assumed that the attenuation is negligibly small, i.e. k = k which 
implies Z s = Z . 

2.4.2 Solutions 



Under the previous assumptions, the standing wave Equation (2.4.2 1 may be 
fulfilled by the real solutions 



with 



<1v(x> t) = (l)o cos(k n x) sin(cu n t) + const, 



2tt 



nA 



(2-4-9) 



(2.4.10a) 



uv = Inni. 



(2.4.10b) 



Here, n denotes a natural number and each of its values stands for a certain mode 
of vibration. Similarly, the value of a> n corresponds to the nth angular resonance 
frequency. 



Using the velocity potential solutions (2.4.9 1, anci in the n g nt of Equations 
(2.3.3 ' anci (2.3.4]!, the acoustic velocity and pressure may be derived as 



v n (x, t) = Vo sin(k n x) sinfcu^t), 



(2.4.11) 



ana 

Pn(x, t) = p COs(k n x) COslcUnt), (2.4.I2) 

respectively, where v and po are the velocity and pressure amplitude related by 

Po = vo Z . (2.-4-13) 
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Chapter \ 
Particle-Image Velocimetry 



The word you've entered 
[velocimetry] isn't in 
the dictionary. 

Merriam- Webster's Online Dictionary 
http://m-w. com 

Despite its high-tech features, Particle Image Velocimetry (PIV) is based 
on a simple idea: recording moving objects at known times and extract- 
ing their displacements by comparing the recordings. Maybe this is the 
most basic understanding one can have of the concept of 'velocity'. 

This chapter is an overview of the principles of PIV and the theory of PIV 
data analysis in general. 

3.1 Historical Overview 

The description of fluid flow can be traced back to very early works of art in 
the way that the artists painted streams, rivers or ocean waves. However, from 
a scientific point of view, maybe the first documented y?cw visualisation was per- 
formed by Leonardo da Vinci in his qualitative studies of cardiac mechanics 
and the vascular system in the 16 th century ([Gharib et al. 2002 1. Figure (3.1 1 



shows Vinci's schematic drawings of the cardiovascular system. 

In 1883 Osborne Reynolds used flow visualisation to demonstrate the lam- 
inar to turbulent transition (Re ynolds|i883 I which enabled him to specify a di- 



mensionless number to describe the flow, later named the Reynolds number. The 
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-|W^* 




I i 

Adopted from {Gharib et al.12002). 

Figure 3.1: Leonardo da Vinci's vision 



flow visualisation employed by Reynolds involved sketching water flow in dif- 
ferent laminar and turbulent phases. Later, at the beginning of the 20 th century, 
Ludwig Prandtl studied the steady and unsteady behaviour of water flow in a 



tunnel by applying mica particles to trace fluid motion at the surface (Prandtl 



The invention of photography provided researchers with the ability to record 
their flow observations with greater reliability and accuracy. It is possible that, 
prior to the scientific application of taking pictures for fluid dynamics research 



purposes, art photographers began capturing images of flows (Grant 1997 1 
Ogle Winston Link is well-known for his photo collection picturing Ameri 



can steam locomotives at the end of their era (Figure (3.2 1). His most renowned 
works were taken at night using elaborate flash bulb systems to capture steam 
particles in air. 

In parallel with progress in photographic technology, namely the increase in 
the spatial resolution of photographs, the innovation of the laser made it possible 
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Figure 3.2: 'Hotshot Eastbound 'by Ogle Winston Link in 1950s. 



to improve the temporal resolution. These improvements paved the way for a 
new family of optically-based flow visualisation measurement systems, with the 
overarching name of laser velocity -meter (velocimetry) . 

Like many other scientific breakthroughs, the first use of laser velocimetry was 
accidental. In their Brownian motion studies using laser scattering, Cummins 



et al. found an unpredicted net shift in the frequency of the incident laser (Cum 
1964)). They created a model calculating the velocity of particles from 



MINS ET AL. 



the Doppler shift and promptly published a paper describing the fundamen- 
tals of Laser Doppler Velocimetry (LDV), also known as Laser Doppler Anemome- 
try (LDA) (JYeh and Cummins||i964>. 



Laser Speckle Velocimetry (LSV), first demonstrated by Barker and Fourney 



( 1977) ), had the novelty of measuring the flow velocity in a less intrusive way than 
that offered by LDV. The technique, based on the measurement of the velocity 
of the visualised fluid-markers, also has the advantage of being whole-field. 

It is difficult to pinpoint the first usage of PIV as this method emerged from 
LSV, borrowing many data acquisition and analysis components from the laser 



velocimetry family. It has been pointed out (Adrian 1991; Westerweel 1993 1 
that when people were performing LSV observations, they were not aware that 
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they were in fact carrying out PIV measurements. An extensive review of the 



development of PIV may be found in Adrian (1986^ 



The theory underpinning the analysis of PIV data is mainly due to Adrian 
( Adrian||i984 ; Adrian and Yao 1985; Adrian|["i988^ ; Keane and Adrian 1990 



1991 1992 1. In its early days, PIV data evaluation was quite time consuming as 
it was based on opto-mechanical techniques. One major advance which has 
eliminated this problem is the use of digital technology in the data process- 



ing, enabling the method to be automated and, almost, real-time ( Willert and 
Gharib||i99i ; Westerweel||i997] ) . 



Nowadays, PIV is a well-established technique. It has been applied to the 
study of flows in many research fields and commercial PIV measurement systems 
are readily available. 

3.2 Pulsed Light Velocimetry 

Consider a system of independently moving objects. An ideal velocity-meter 
which can measure both the magnitude and direction of the velocities of the 
objects should follow the rationale below: 

data acquisition. Marking the moving objects in space and time. 

Rule 3.1. Spatial marking is necessary when the actual objects are not vis- 
ible. It can be achieved by adding some easy-to-visualise but faithful-to- 
follow external objects to the system. 

Rule 3.11. The spatial marking must not be intrusive to the system. In 
addition, the markers must not interact with each other. 

Rule 3-iii. Temporal marking is defined as capturing the spatial state of the 
system at a given time. The marking time should be short enough so that 
the moving objects can be considered as 'frozen' within it. 

data analysis. Evaluating the recorded markings. Based on the funda- 
mental definition of velocity and depending on the characteristics of the 
marking data, an evaluation process can be applied to yield the velocity 
vectors corresponding to the observed objects. 
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This abstract velocity-meter can be applied to fluid flows almost without los- 
ing its generalisation. Since the fluid elements under investigation are not visible 
normally, some fluid markers, also called tracer particles or seeding particles, which 



fulfill rules <3.i) and ( 3 .ii > are added to it. Based on rule < 3 .iii > , the tracer parti- 
cles are illuminated at least twice for a short enough time and their images are 
recorded for subsequent processing. The analysis process converts the images 
to a grid of velocity vectors or a vector-map. As the temporal marking part of 
the measurement is usually carried out by short pulses of light, the technique is 
referred to as Pulsed Light Velocimetry (PLV). 

3.2.1 The PLV Class 



Adrian (1991I classified PLV according to the attributes of the applied fluid 



markers and the phases of the PLV images (Figure (3.3!). In the first step, the 
PLV class is organised according to the fluid marker type. Molecular markers, 
which are activated by laser beams to visualise the flow, are used only in the 
study of certain flows whereas particulate markers have a more general usage. It 
is particulate markers which are of particular interest in this study. 

Depending on how the PLV data collection system is configured, for par- 
ticulate markers the acquired data falls into one of three phases, each of which 
requires its own specific method of analysis. These three phases are explained in 
the next section. 

3.2.2 PLV Data Phases 

Raw PLV data comprises an image which records the positions of the fluid mark- 
ers, referred to as particle-images, at an arbitrary time. The phase of a PLV image 
is specified by two dimensionless numbers. We review these numbers here and 
leave their analytical definitions for later. 

source density !Ns- A number that determines whether the particle- 
images are distinct individuals or whether they overlap each other in the 
image plane and form a speckle. 

image density Ji \. In processed PLV data, which has a grid structure with 
velocity vectors at each grid point, a vector is representative of a specific re- 
gion around it known as an interrogation area. Image density is the number 
of particle-images within an interrogation area. 
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Molecular Markers 




Photochr 



f . Pulsed Light Velocimetry (PLV) \<f ~'^\ Fluorescent ] 

{Particulate Markers) 




Particle-Images (PIV) 



Speckle Patterns (LSV) 
X 5 > 1 




Low Image Density PLV (PTV) 
N/< 1 



High Image Density (PIV) 
W/> 1 



Figure 3.3: Pulsed Light Velocimetry class as grouped by Adrian (ippi). 



Different values of source and image densities result in three possible phases 



for PLV data (Figure (3.4!) 



3.2.2.1 Laser-Speckle Phase: N s > 1, 3^ > 1 

When the source density DSf s is too high, which obviously requires the image 
density !Ni to be high too, the particle-images cannot be found in the single 
form. Instead, the light scattering from the tracer particles overlaps and creates 



a pattern known as laser speckle (Figure (3.4c)). The image processing method 



require for this phase determines the displacement of speckles between successive 
exposures, a technique referred to as Laser Speckle Velocimetry (LSV). 

3.2.2.2 Low-Density Particle-Image Phase: N s < 1, Ni <C 1 

In this case the source density !N s is so low that the distances between particle- 
images are much more than their dimensions. This leaves most of the interroga- 
tion areas empty. The ones which do contain some scattering particles have so 
few of them that it becomes easy to link the particle-images from two successive 
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a) PTV:Ji s <. 1.N/-C 1. b) P/K-K 5 < 1,K 7 » 1. c) Z5V/X 5 > 1,X 7 » 1. 

77?<? circles display interrogation areas. 

Figure 3.4: PLV data phases. 



exposures; i.e. it is practical to 'track' each individual particle (Figure (3.4a)). The 
related image analysis method for this phase is called Particle Tracking Velocime- 
try (PTV). Because of the lack of spatial data, it is not possible to get a regular 
multi-field grid of velocity vectors from this particle-tracking' phase. 

3.2.2.3 High-Density Particle-Image Phase: 3\f s < 1, Ni > 1 

This phase falls between the laser-speckle and particle-tracking phases. That is, 
the source density IN" s is low enough to give separate particle-images, but there 
are sufficient scattering particles contained within each interrogation area for 
velocity vectors to be extracted (Figure ( 3.4b! ). The high image density 3Nf 1 in- 
creases the ambiguity in matching the scattering particles in successive exposures, 
making particle-tracking style processing inefficient. Instead, the displacement 
of particle-images is calculated 'in group' using statistical techniques. This is the 
essence of Particle Image Velocimetry (PIV) image evaluation. 



3.2.3 PIV Principles 

The remainder of this chapter is focused on the high-density particle-image, or 
PIV, phase. The principles of PIV are, in general, the same as those described in 



the implementation of the abstract velocity-meter at the beginning of §3.2 The 
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forthcoming sections, as outlined in Figure (3.5!, discuss each element of the PIV 
method in detail. 



[Flo 




I Data Acquisition \- ^ _ __ 



Post-Processing 



A V 



Visualisation J 



* -[ Particle-Imaging 



■-G 



Re con 



Figure 3.$: Principles of PIV. 



3.3 Visualisation 

3.3.1 Light Sources 

While white light sources have the advantages of being low-cost and not re- 
stricted by laser safety rules, they are not the best choice as they are less colli- 
mated than monochromatic sources and recording mediums are not optimised 
for them. 

Lasers have been used widely in PIV set-ups since the early days of the tech- 
nique. Apart from being free of the disadvantages of white light sources, the 
variety of laser types makes them suitable for specific purposes. 

cw. Continuous Wave (CW) lasers have an output beam which can be oper- 
ated continuously. The power is adjustable and, by applying a mechanical 
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shutter, it is possible to generate short pulses from the output. However, 
the output power may not be high enough for high-speed flow visualisation 
or short-exposure PIV recording. 

pulsed. Unlike continuous wave lasers, pulsed lasers have a discrete out- 
put. Pulsed lasers can be one or more of the following types: 

q-switched. The output of a Q-switched laser is a short-duration 
high-energy pulse which is generated using a shutter-like device. I 
The device releases the output only when the built-up laser energy 
reaches a specific limit. 

single-pulsed. Also known as the normal mode, the output of a 
single pulsed laser consists of only one pulse which is of the order of 
0.1 to 1 ms. 

repetitive-pulsed. Also known as scanning lasers, the repetitive- 
pulsed lasers have an output of the order of 10 4 pulses per second. 

mode-locked. Although laser beams are monochromatic, in the laser 
cavity there is a range of modes close in frequency to each other; this 
can affect the characteristics of the output. The locking technique 
synchronises the phases of all the modes and, consequently, generates 
an output of regularly-spaced pulses. The duration of each pulse is 
around 10 15 to 10 12 s. 



Table (3.3.1! lists some common lasers employed in PIV and describes their 
properties. It can be seen that Nd:Yag lasers have features which are very appro- 
priate for PIV measurement. 

3.3.2 Light Sheet 

The minimum optical configuration needed to create a sheet from a beam of 
light is a cylindrical lens. When the light beam is narrow and collimated, a 
single cylindrical lens can be effective in generating a sufficiently thin light sheet. 
However, in most cases it is preferable to add a few more optical elements into the 
configuration to allow control of the thickness, height and intensity distribution 
of the sheet. In this section, various ways of generating light sheets, along with a 
formulation of their intensity distribution, are presented. 



'The Q" refers to the measure of the sharpness of the resonance peak. 
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Laser 


Name 


Wavelength 


Power 


cw 




Pulsed 








urn 






QS* 


MLt 


SP* 


RP S 


He-Ne 


Helium-Neon 


633 


10~ 3 - 10~ 2 w 


/ 










Cu 


Copper-Vapor 


510,578 


1 -30 W 








/ 




Ar+ 


Argon-Ion 


488,514 


60 - 100 W 


/ 


/ 




/ 




Cr+ 3 


Ruby 


694 






/ 




/ 




Nd:YAG 


NeodynrYAG 11 


532 


> 0.4 J per pulse 


/ 




/ 




/ 



Q-Switch. 
1 Mode-Locked. 
* Single-Pulsed. 

Repetitive-Pulsed. 
™ Yttrium-Aluminium-Garnet. 



Source: Young (ip86) 



Table 3.1: General properties of some lasers used in PIV. 



3.3.2.1 Configurations 



Figure \^.6\ shows some schematic optical configurations for creating light sheets, 



as suggested by Raffel et al. (1998 ). 

One way to manipulate the light sheet emerging from a cylindrical lens is to 
add two more lenses of the same type in order to vary the height and thickness 



(Figure (3.6a!). Because of safety reasons, which will be explained shortly, the fo- 



cal point of this configuration is chosen to be virtual. This restricts the structure 
spatially and consequently puts a lower limit on the thickness of the final light 
sheet. 

Sandwiching the sheet-generating cylindrical lens between two spherical lenses 
(one diverging and one converging) provides the option of more dynamic vari- 



ation of the height and thickness of the sheet (Figure (3.6b)); nonetheless, these 



parameters cannot be controlled independently of each other. The configuration 
has the advantage of being more economical by using spherical lenses as opposed 
to cylindrical. 

Both of the previous configurations contain a lens with a negative focal length. 
This ensures there is no real focal point in the set-up which prevents any high- 
energy concentration point in the system. When using a high-energy laser as 
a light source, focal points are generally to be avoided as they can burn small 
particles like dust floating around the point or, in a less likely case, ionise the 
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Side View 



Top View 




-50 mm +200 mm +500 mm 



A 



Azz 



a) —Cylindrical, + Cylindrical, + Cylindrical 



View 



Top View 




b) —Spherical, -\-Cylindrical, -\-Spherical 



View 



Top View 




c) -\-Cylindrical, + Cylindrical, + Cylindrical 

Shifting the lenses respect to each other adjusts the height and/or thickness of the light sheet. Adopted 
from Raff el et al { ipp8 ). 



Figure 3.6: Some light sheet generation configurations. 



air particles. In cases where versatility is crucial, a configuration such as the one 



depicted in Figure (3.6c) can be used in conjugation with appropriate safety pre- 
cautions. 
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3.3.2.2 Intensity Profile 

If the generated light sheet is assumed to lie in the x-y plane, and the illuminated 
area of interest is small enough to make the light intensity distribution constant 
in the plane, the intensity profile will be a function only of the depth z. The 
spatial variation in intensity of light emitted from a point source is Gaussian in 
nature; the same is true of light from a line source. Hence, the intensity profile 
of a light sheet can be expressed as a Gaussian function of z: 

J(z)=Joe- 8iz - Zo]2/A ^ 2 , (3.3.1) 

with J > zo and Azl being the maximum intensity, position of the centre and 
thickness of the light sheet respectively. 

Some modern optical set-ups offer an intensity profile with a rectangular box 
form given by 

J Jo if |z-z | < Az L /2 

elsewhere 



3.4 Particle-Imaging 



Figure (3.7 1 illustrates the optical details of the particle-imaging system in a PIV 
set-up. The camera's convex spherical lens is assumed to be thin and aberration- 
free 2 with a focal length of Yf and diameter of D L . The average distance of the 



light sheet described in §3.3.2 to the camera lens is d , while the average distance 
of the camera lens to the image plane is d[. The finite thickness of the light sheet 
Azl causes the visualised seeding particles to lie at different depths from the lens. 
Therefore, their images are magnified by a factor which depends on the depth z. 
Assuming the origin of the axes is chosen to be at the centre of the light sheet: 

, , r -. d; Azl Azl / s 

M(z=— — , -— ^<z<+— \ (3.4.1) 

d — z 2 2 

where d Q and d; are related to each other by 

111 , . 

X- + T = -- (3-4-2) 

d di yf 

2 ^3-4-5 
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Image Plane 




Figure 3.7: Particle imaging system. 



Throughout, it is assumed the light sheet thickness is much smaller than the 
distance of particles to their images, i.e. Azl <C (d D + di). This implies that 
perspective projection can be ignored; that is, the tracer particles can be mapped 
from the light sheet volume to the x-y plane, referred as the object plane. The 
mapping for any function f( x) is given by its mean over the intensity across the 
light thickness: 

1 



r + Az L /2 



<f(x)> 



J(z) 



% 



J(z)f(x)dz. 



(3-4-3) 



-Az L /2 



The lack of perspective projection means that the variation of the magnification 
M(z) across the light sheet thickness is negligible and can be written as 



M = M(0) 



do' 



(3-4-4) 



The purpose of the particle-imaging system is to map the positions of the 
tracers in the object space (the light sheet) to the recording media image plane. 
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This transformation involves some loss of data originating from both the imaging 
and recording processes. The aim is to anticipate the output of the imaging 
device considering the data loss error. To facilitate this, diffraction theory is 
used to describe the formation of particle-images and then their diameters are 
determined as a function of the imaging system parameters. 

3.4.1 Scattered Light Field 

In this section, the analytical structure of the light field is discussed. Through- 
out, it is assumed that the light source used for illuminating seeding particles is 
monochromatic and coherent. In common with the particle-image expression, 
the corresponding tracer particle is sometimes referred as a particle-source. 

Let us consider the scattered light wave from the particle-source, a{x), to 
have an angular frequency cb and complex field amplitude 2L[x). The field can 
be described by the following equations (Goodman||i968, §3.2) 



a x 



m[Mx)e-^], (3.4.5a) 



Mt) = Mt)e- x ^ {lt] ; (3.4.5b) 

where S%{ x) and cp( x) are a non-imaginary amplitude and phase respectively 3 
In this context, it is the time-invariant part of the wave, the complex field ampli- 
tude, which is of interest and which satisfies the Helmholtz equation 

(V 2 + k 2 )^(x) = 0, (3.4.6) 

where k is the wave number, ordinarily related to the wavelength A and group 
speed c by 

k = 27t/A = 27tcb/c. (3.4.7) 



3.4.2 Diffraction Pattern 

In spite of the geometrical optics prediction, typical particle-images on a record- 
ing medium plane are not regions with sharp borders that are well-separated from 



3 To distinguish between light field and sound field, the parameters for light waves are marked 
by check signs ( " ). 
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the background; instead, they fade out to it in the form of fringes. This can be 
explained by Huygens's principle founded on the wave nature of light. 

Particle-images can be thought of as a colony of diffraction patterns, pro- 
duced by light waves scattered from the tracer particles, which are incident on 
the imaging system aperture. Depending on the distance of the tracer particle 
from the aperture d , the distance of the aperture to the image plane d;, the size 
of the aperture (camera lens diameter) D L , and the wavelength of the scattered 
light A, the produced diffraction pattern can be grouped into different regimes. 
The regimes can be specified by a dimensionless quantity known as the Fresnel 
number: 

^ = 7^-. (34.8) 

4 A d 

Generally, in a PIV experiment, the Fresnel number is much smaller than unity, 

N F <1. (3.4-9) 

This means that the object distance d D is much greater than the aperture size, 
implying that both the incident and diffracted light waves are effectively planar. 
Such a regime is known as Fraunhofer or far-field diffraction. 

The diffraction pattern of a particle-image can be treated as a two-dimensional 
distribution of irradiance 4 in the image plane. In the case of a monochromatic 
field, the irradiance distribution is simply the square of magnitude of the corre- 
sponding complex amplitude, 

o(l) = \Mt)\ 2 . (3-4.10) 

The form of the irradiance distribution can be determined by the diffraction 
pattern regime and the aperture properties. For a common PIV set-up, the dif- 
fraction pattern, as mentioned earlier, is of the Fraunhofer type and the camera 
lens is a circular aperture. Bearing these in mind, and using the Kirchhoff-Fresnel 
formula, which is a precise mathematical extension of Huygens's principle, it can 
be shown that the irradiance distribution of the particle-images takes the form of 
an Airy pattern ( Hardy and Perrin|i932 , §60). This pattern is expressed by the 



4 Sometimes intensity is used as a synonym of irradiance or illuminance . As discussed in 



Fowles (1975 §7.2), the word 'intensity' is not technically correct in this context. Hence, unlike 
the other resources, usage of the expression image intensity field in this work is avoided and 
substituted by irradiance distribution. 
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Airy function which in the polar co-ordinate system of the particle-image is: 



1,-1 i / -ft( f / f o) 'i M 
J l r J = J o ( —^ — 1 , with 



(3.4.11a) 



D, 



2f 



(3.4.11b) 



TO 



Ad; _ 

7tD L ' 



(3.4.11c) 



and j?i (f) denotes the first order Bessel function of the first kind. Equation 



(3.4.11a! describes a circularly symmetrical diffraction pattern consisting of a 
bright central disk surrounded by a set of dark-light concentric annular bands. 
The set starts with a dark band and the illuminance of the light bands diminishes 



rapidly. The bright central area is known as the Airy disk (Figure (3.8)) 



logj(f) 




Figure 3.8: The Airy pattern irradiance distribution in particle -image plane. 



30 



3-4- Particle-Imaging 

3.4.3 Impulse Response 

An imaging device can be considered to be a linear system with both an input 
and output. To measure the 'quality' of the system, a standard input is needed 
as a basis for comparison. The standard input normally used is a sudden surge, 
or impulse, with infinite amplitude and zero width. This surge, better known as 
the Dirac 6-function, is referred to as the unit impulse. The output of a system 
when the input is a unit impulse is the impulse response or point-spread function 
of that system. 

The impulse response may be treated as a 'mapper' which transforms the 
particle-source field amplitude in the object space S^[x) to the particle-image 
field amplitude in the image plane Ai{x). This mapping is given by the super- 
position integral 

* r+00 

^i(x)= 3iAx\ x)S^[x)Ax, (3-4- 12 ) 

j j —00 

where the comp lex impulse response is denoted by !Ka( x; ~x). The integral in 



Equation (3.4.12 1 is said to be the convolution of the functions Ji A [ x ) and S^{ x 



and in a simplified form is shown as 

Ai(%) = X*[t) © -3oU), (3-4-I3) 

where the sign © denotes the convolution operator. Broadly speaking, the more 
the input J4„ and the output S\ look like each other, the higher the quality of the 
imaging system. 

In PIV theory, it is assumed that the particle-sources are geometrical points 
described by the Dirac 6-function, i.e. the input is a unit impulse. As a con- 
sequence, the output is simply the input impulse response of the system. That 



is, the impulse response for a PIV system is given by the Airy pattern (3.4.11a! 
which in its radial form in the coordinate system of the particle-image is given 
by 



r,M2 W) 



|tt*W = -^. (3-4-14) 

Jo 

This is an important property of the system which is not often explicitly men- 
tioned in the literature. 



Practically, when used in calculations, equation (3.4.14! is substituted by a 

3i 



Particle-Image Velocimetry 



Gaussian approximation ( Westerweel|i993 §A.4) given by 

|^(f)| 2 = e- (f/2f ° )2 . 



(3-4-I5) 



3.4.4 Source-Image Mapping 

It is convenient to break up the mapping of particle-sources from the visualisa- 
tion volume to the recording plane into a few steps. 

The first mapping is from the object space to the image plane using the geo- 
metrical optics theory. Consider, for the moment, that the diffraction effect in an 
imaging system with a converging lens is negligible. Under this condition, the 
prediction provided by geometrical optics theory is predominant. That is, the 
particle-image is the exact reproduction of the particle-source, but, magnified 
and inverted; 

x p = -Mx, y p = -Mu, (3-4-i6) 

where x p denotes the position of the particle-image as a geometrical point. 



A similar relationship exists between the complex field amplitudes (Goodman 

1 



1968 §5.3) 



A: 



(G) 



M 



-3oU 



(3-4-I7) 



where SV indicates the geometrical prediction. 

Finally, to obtain the actual field amplitude of particle-images, the diffrac- 
tion effect must be incorporated. This is a mapping from the geometrical point 



form of the particle-image to its diffraction pattern. Goodman (1968) gives the 



final image output as a convolution of the aperture impulse response and the 
geometrical prediction. The convolution in our notation can be re-written as 



S\{x 



■> + oo 



MJflA 1 



dx p , 



(3.4.18) 



where as illustrated by Figure (3.8 1 



Substituting Equation (3.4.17) into the convolution expression gives the source - 
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image mapping in the following compact form 

1 



M* 



M 



CK; 



■3o(x 



(3-4-19) 



The convolution (3.4.19 1 expresses the final form of the complex amplitude 
field in the image plane, however, in practice the quantity of interest is the ir- 
radiance 5(x). For any given position x in the image plane, the square of the 
magnitude of the complex amplitude field generated by two particle-images po- 



sitioned at x n and x ' is 



w« 



p 



^(x 



pj •* '"SAX 






(G) (? D U (G) *( 



Xp)dx p dXp. 

(3.4.20) 
Assuming the imaging system is set to be in the PIV phase, the source density 
IN" s is very low. Hence, there is no overlapping between neighbouring particle- 
images. This implies that, at each position in the image plane, the irradiance 
value is generated by only one particle-image or 



4 G) U p )^ G, *(? p ) 



14 G) (x p )| 2 6(x p -x. 



This greatly simplifies Equation (3.4.20) to 

|^U)| 2 H^(xV©l^ G) U)l 2 , 

which can be written alternatively as 

1 



J(x 



M 2 



l^(xr©j(x). 



(3.4.21) 



(3.4.22) 



(3-4-23) 



In other words, in the PIV phase, the impulse response of the light intensity 
distribution is simply the square of the magnitude of the imaging system impulse 
response. 

3.4.5 Resolution Limit of Camera Lens 

The Fraunhofer diffraction imposes a resolution limit on the imaging aperture. 
Optical systems with this limit are said to be diffraction-limited or aberration-free. 
According to the Rayleigh criterion, the minimum distance between two particle- 
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images at which they can still be 'resolved' by the camera lens is the radius of 
the corresponding maximum Airy disk (Hardy and Perrin 1932, §61). This 



distance is known as the theoretical resolution limit. 

The Airy disk radius, ta, can be obtained by calculating the first zero of 
the Bessel function J?i(f/f ) occurring for f/f = 3.8217. Applying Equation 



(3.4.11c!, the radius is given by 



f A = 1.22A 






(3-4-24) 



It is convenient to express this quantity as a pure function of the system pa- 
rameters. Combining Equations (3.4.1! and ( 3.4.2] ), Equation (3.4.24! can be 
alternatively written as 



A L = f A = 1.22A (M+ 1) f # ; 



(3-4-^5) 



where A L denotes the resolution limit of the camera lens and f # is the f-number 
of the aperture lens defined as 

c# = JL 



f 



(3.4.26) 



3.4.6 Particle-Image Size 



The convolution equation (3.4.19 ) together with the system resolving limitations 



can be used to calculate the particle-image dimensions. Adrian and Yao (1985) 



used a Gaussian approximation for the convolution components — the Airy irra- 
diance distribution and the image predicted by geometrical optics — and obtained 
the following estimated value for the particle-image radius 



M 2 r p 2 + ^A L 2 + ^A R 2X 



where A R notates the resolution limit of the recording medium. 



3.5 Spatial Correlation Theory 

The concepts presented in this section are based on appendix |A] 



(3-4-27) 
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3.5.1 Source and Image Densities 



Earlier, in §3.2.2 the concepts of source and image densities were introduced. In 



this section, more accurate definitions for these quantities, as given by Adrian 



(1984 1991I, are presented. 

Consider a cylindrical volume in the particle-source space which is created 
by the intersection of the illuminating light sheet (of thickness Azl) with a circle 
whose radius is f p /M; i.e., the circle that is produced when a particle-image with 



a radius given by Equation (3.4.27 1 is back-projected on to the object space. This 



cylindrical volume is said to be a resolution cell expressed by 

VR = AZL7t (M) 2 - M 

The mean number of particle-sources in a resolution cell is the source density 
and is given by 

Ws = CV R , (3-5-2) 

where C is the mean number of particle-sources per unit volume. 

From the definition of the source density, it can be immediately concluded 
that the condition for particle-images to be thoroughly separated from each other 
is that the corresponding particle-sources are laid uniquely within resolution 
cells. That is, if more than one particle-source shares the same resolution cell 



their associated particle-images overlap. Applying Equation (3.5.2), the probabil- 
ity of exactly k particles being in a volume Vr is given by the Poisson distribution 
( |A.2.i6| ) as 



For N s > 1, the probability p- k>1 . yi becomes significant which means forma- 
tion of speckle patterns, whilst for J\f s *C 1 one has 

, :N s k 

implying the phase is particle-image. 

By analogy with the resolution cell definition, an interrogation cell can be 
defined as 

V I = Az L 7t(^) 2 , (3.5.4) 
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which describes a cylindrical volume created by the intersection of the light sheet 
with a circle whose radius is given by back-projecting an interrogation spot with 
radius Di in the image plane. The mean number of tracer particles in an interro- 
gation cell is the image density given by 

Ni = CVj. (3.5.5) 

In the case that Ni » 1, the probability p. , N is high, meaning that there 
are many particle-images in each interrogation spot; i.e. the phase is high-density 
particle-image. On the other hand when Ni < 1, the probability P- k<1 .- N is 
more significant implying the phase is low-density particle-image. 

3.5.2 Displacement Field 

Assume that the flow under study is described by the velocity field v" I x"(t) 1 . 
If an ideal tracer particle at time t t is located at position x"(t) with velocity 
v~p I x"(t) 1 , based on Rule (3J), at time tj the particle displacement field is given 

by 

ft; 

(x(t)j dt. (3.5.6) 



V P 



d?-)i = ^(x;tt,tj) = 

The displacement field is an 'averaging' over time, which means it does not 
give any information on the motion of the flow within the averaging time. Thus, 
there always exists some error e where 

kx^-i- v* fx"(t)J (tj-tOl^e, ti^t^tj. (3.5.7) 

Needless to say, a PIV experiment only measures the velocity of the flow 
where the tracer particles are positioned. That is, since the particles are located 
randomly, the flow can only be probed by random sampling. 

3.5.3 Particle-Source Pattern 

Particle-sources can be treated as geometrical points randomly distributed in the 
object space. Assume there are N p tracer particles in the flow under investigation 
and at any time t the particle a is positioned at x" a (t). The scattered light 
intensity, at any point x"(t) in the object plane, may be described by the Dirac 
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6-function as 

, J{x)=J^J(z 0C )8{x-x a ,y-y a ). (3.5.8) 

a=l 

This equation can be given in a compact form by 

(x(t)) =J (g (x(t))) (3.5.9) 



3(z 

where 



G (x(t)) =2^8 (?(t)-? a (t)} . (3.5.10) 



cx=l 



The dimensionless quantity G I x"(t) I is known as the particle-source pattern. 

For each element t t of a given time set {ti,t 2 ,. . . , t n }, the particle-source 
pattern has a value G I x[t\) I which for brevity may be denoted as 

G t (x) = G i (x(t)) =G [x(tOJ . (3.5-11) 

As time passes, from all of the possible patterns (the ensemble of patterns), a 
'history' is generated which describes a stochastic process? We describe the sto- 
chastic process of the particle-source patterns as 

{G i[x),G 2 {x ),..., G^x)}, 

or in a compact form as {G( x")}, where G( x) is said to be a random variable 6 or 
a realisation. 

Let T^.jt indicate the displacement between two points [x, t t ) and [x 1 , tj) 
by 

t^)x=x'{t ] )-x[t i ], (3-5- 12 ) 

The process {G(x")} may be considered as a chain of events in which each indi- 
vidual random variable Gj( x + TVjjt) at ( x'+TV.^tj), as an output, is produced 
by another random variable Gt( It) at ( x", t^), as the input. In order to link these 
input and output signals, an impulse response ^KGjGi(?V ; ji) can be specified and, 



6 



5 £A.4 



Definition (A.vl 
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consequently, each random variable can be written as the following convolution 



Gj( s Tt-ji.) 



^GjGj^^jO © Gi(7^;ji] 



(3-5-I3) 



Assuming the displacement field to be local, i.e. cC^i = d ^, the impulse re- 
sponse can be expressed as the shift of the input by the local displacement: 



^GjGil * x;ji) — 5( * x;ji — ^jij! 



(3-5-I4) 



At any point x", the data of velocimetry interest to extract from the process 
{G( x)} is the displacement ~tt,)i- One way to do this is to find out how 'corre- 
lated' the random variables Gt( x) and G 3 -( x + "?V ; ji) are by calculating the co- 



variance matrix 7 corresponding to the process {G( x)}. Using Equation (A.4.18 



the covariance matrix elements for the random variables are given by the follow- 
ing cross-covariance 8 



X 



GjGil ■> x;jij 



cov{G i (x + ),G j (x > + 7^- n )} 



= £{ (g^x) - £{G}(x)}) (g^x + 7f,n) " E{G t (x + ?*;,i)}) }, 

(3-5-I5) 
where cov{ , } stands for the covariance operator 9 and the expectation value IO is 



denoted by £{}. Substituting Equation (3.5.13! into the cross-covariance relation 
yields 



%. 



GjGjl ■> x;jij 



cov{Gi(x"), 5C G jG t (x"+ 7 



x;jij 



Gj(x"+ '/?.«)} 



^GjG l ("j'^;ji) ©COv{G l (5c"),G j (5c* + l*V;ji)} 



(3.5.16) 



= ?t" GjGi (lV;ji) ©^ GtGt C?*^;ii), 

Defining the auto-covariance a values as the diagonal members of the covariance 



76A. 

10 tX. 



4-3| 

4-2-1 
4.1.2 

4.1.3 
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matrix denoted by 



%. 



G^ •> x;ij 



*L 



GjGil ■> x;iij 



the cross-covariance can be alternatively given by 

^G j G l (^ ; ij)={^G j G i ©^G t }(^;i)- (3-5-I7) 

That is, the non-diagonal members of the covariance can be generated by know- 



ing the diagonal members. Using Equation (3.5.17), it can be shown (Adrian 
1988^ that the covariance matrix elements are of the form 



%. GjGil s 5?;ji) 



Ch{~s^.»- d n ) 



(3.5.18) 



Each element of the covariance matrix given for a particle-source pattern de- 
scribes a point in the displacement-probability space; i.e. ^ggJ^Vjji) repre- 
sents the likelihood that the displacement of particles from the realisation Gt( x") 
to the realisation Gj(x" + "?V;ji) is equal to to d ^. The set of these points is 
referred to as the correlation plane. 



The covariance matrix in Equation (3.5.15! describes a correlation plane with 



one infinite sharp peak and zero elsewhere. In the velocimetry of a particle-source 
pattern process by an imaging system, the aim is to 'measure' the correlation 
plane of the process. That is, to find the correlation plane corresponding to 
particle-images. Due to the limiting errors in the imaging system, discussed in 



§3.4.3: the measured correlation plane contains some noise and has a wider and 



shorter peak. In the next section the formulation of the measured correlation 
plane is presented. 



3.5.4 Particle- Image Pattern 



Earlier, in §3.4.4 the impulse response of light intensity was described. The 



impulse response may be employed to derive irradiance distribution from light 



intensity distribution. Combining Equations (3.4.23) and (3.5.9), the following 
convolution exists 



5(x)=5o|^(x)| : 



1 



<G(x)> 



3( z ; 



(3-5-I9) 
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By analogy with the particle-source pattern definition, the second part of the 
convolution is named the particle-image pattern and is given by 



G(x) = (G(x,-Q,z)) a(2) , 
or, alternatively, in the explicit form of 



(3.5.20a) 



N, 



G(x) = 3- ^5(z a )6(x - xj 
Jo 



(3.5.20b) 



<x=l 



Let d i) and 5 ^ be, correspondingly, the equivalents of d ^ and l*V;ij in 
the image plane given by 

d n = M[9..7 n ,i}.7 n ), (3-5-2i) 

5 *;ji = M (£-^V;ji.D-^V;ji)> (3-5-22) 

where £ and i) are the unit vectors and . denotes the scalar product operator. 
Knowing the source-particle pattern form, the covariance matrix of the particle- 
image pattern process can be obtained ( |Adrian|i988£ |. The cross-covariance of 
two realisations G\ and Gj is of the form of 



^G 



GjGjV S t;jV 



O 



c 

M 2 



H 5 x;]i — dji), 



(3-5-23) 



where C o is the number of particle-sources per unit area which remain within 
the light sheet expressed as 



Co = CM 2 Az L F o (£.^0- 
Here Fq(z. d ^) is the out-of-plane loss-of-correlation given by 



(3-5-24) 



Y [z.d n ) 



3{z + z.d 



n ■ 



Hz) 



mz)) 



(3-5-25) 



0(z) 



Comparing Equations (3.5.18 1 and (3.5.23) indicates that the 'measured' ^.cg^ 



or $t q.q , has the same structure as the original ^. G . G .; but mapped within the 
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imaging system and altered by a constant. The constant involves number of 
source particle-sources in the illuminated volume and an error originating from 
particle-sources losing their pairs between the successive realisations. 

As the final step in calculating the spatial correlation properties of PIV data, 
the cross-covariance of the irradiance can be derived using the known cross- 



covariance of the particle-image pattern. This derivation as given by Adrian 



(1988^1 may be written down as 



with 



anc 



fio 



p + 00 



|UUx)| 2 dx, 



"JfU? 



{|JG 



l^l 2 }U? ;n ) 



n' 



% 



(3.5.26) 



(3-5-27) 



(3.5.28) 



as the normalised self-correlation peak. Substituting Equation (3.5.23 1 into (3.5.26) 
yields 

^.5 J if l (^ : ji) = Co5g^F M (Tit. Jt - i n ). (3.5.29) 

In its new form, the correlation is an indirect function of the complex amplitude 
impulse response. This implies the correlation peak width is proportional to the 
dimensions of the particle-images. 



3.5.5 Interrogation 

The PIV correlation theory presented so far describes the displacement of flow 
within two realisations without any limited spatial boundary. In practice, it is 
desirable to measure the displacements inside the interrogation cells. To do this, 
the irradiance field in each realisation is limited to a region of interest defined by 
a window function. 



Like the light sheet intensity profile, described in §3.3.2.2 the typical window 
functions are of the Gaussian or rectangular box form given by 



•win 



W e 



— 8 X — Xi 



iVA^-SH-ftjVAft 2 



(3-5-30) 
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or 



<W if |x-XiKAxi/2 A |^ — -Qil < AfJi/2 

■^(x)H n 1 U ' (3-5-31) 

elsewhere 

where (xi,tji) and (Axi,Ayi) are, respectively, the centre and dimensions of the 
interrogation area. For consistency, throughout, the interrogation area is denoted 
as Si. 

The interrogated irradiance field is defined as 

i(?) = ^(?)5(?). (3.5.32) 

Given 

I t (?) = I (?(t0) , (3-5-33) 

the cross-covariance matrix for I t ( x ) and Ij( x ) over a shift ?V ; ji is given by 

1Li } i l (?t;n) =cov{I i (?) ) Ij(?+T^.. i )}, (3-5-34) 



which may be calculated by a similar procedure presented in §3.5.4, Introduc- 
ing the window function changes the averaging of irradiance distribution, conse- 
quently, some new terms appear in the outcome cross-covariance. It is convenient 
to present the cross-covariance in terms of its ensemble mean and fluctuations by 

£n(?) = (*. n (?)) + £n(?), (3-5-35) 



where, for the sake of brevity, the realisation information indices are omitted. In 
the next section, the ensemble mean covariance is studied in detail. 



3.5.6 Analysis Modes 

Consider a frame is a recording medium at a specific time, and an exposure is 
registering a realisation of the image plane on a frame. PIV data acquisition 
has the option of recording images with a combination of number of frames 
and exposures. Depending on this combination, the ensemble mean covariance 



appeared in Equation (3.5.35) has different forms as presented in this section. 
42 



3- 5- Spatial Correlation Theory 



3.5.6.1 Single-Frame, Single-Exposure 

While no velocimetry can be carried out in this analysis mode, it includes some 
components which can be used in the other analysis modes. As there is only one 
frame, the ensemble mean covariance is, actually, auto-covariance which in its 
normalised form is known as auto -correlation. a In this special case, since there 
is only one realisation, the auto-correlation is a self-correlation. Following the 
Adrian[ (1988^), tne auto-covariance % l can be decomposed to the 



notation in 
following 



%. 



%. 



d 



+ £1 



+ £1 



(3-5-36) 



where % c stands for the mean background correlation, %. F is the fluctuation of 
the background noise and %. P shows the self-correlation peak, also known as the 



pedestal (Figure (3.9)). As expected, the pedestal is located at the origin of the 
correlation plane. 




Figure j.p: Typical correlation plane for single-frame single-exposure mode. 



^A.4.1.3 
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3.5.6.2 Multiple Frames - Single Exposure 

The ensemble mean covariance in this mode is cross-covariance or, when it is 



normalised, cross-correlation.^ Keane and Adrian (1992) gives the double-frames 
single-exposure ensemble mean cross-covariance as 



£ 



in 



£cU) + £fU) + £dU 



(3-5-37) 



Comparing this covariance to the auto-covariance Equation (3.5.36 1, the pedestal 




Figure 3.10: Typical correlation plane for multiple-frames single-exposure mode. 

is replaced by the displacement-correlation peak 3^ D which contains the informa- 
tion of the mean displacement of the particle images within the interrogation 
area (Figure ( 3.10) )). The displacement-correlation peak, in its precise form, is 
given by 

*d( S? ; ,t) = <?oi( sV ;ji ) % hi F M (?^ t - d n ). (3.5.38) 

This relationship is comparable to the cross-covariance of the interrogated irra- 



diance fields given in Equation (3.5.29 1. The only difference is that the latter 

I 3^A.4.2.I 
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is described by the number of particle-images per unit area which stay within 
the light sheet and the interrogation areas between successive realisations. This 
quantity, denoted by C o\, is referred as the effective particle-image density and is 
given by 

coi(T^ t ) = ^iF I (T^, l )Fo(£.7, t ) J (3.5-39) 

with 

W ?.„) = ^-{Wj ® WiX ??.„), (3-5.40) 

which is known as the in-plane loss-of-correlation. 



3.5.6.3 Single Frame - Multiple Exposures 

When, due to hardware limitations, the image recording process is not fast 
enough to capture each realisation in a separate frame, the resultant data can 
be analysed in a way which is a combination of the previous modes. The out- 
come irradiance field, denoted by I z, is a super-position of the individual reali- 
sations. As there is only one frame, the displacement can be extracted from an 
auto-covariance. In the case of a single-frame double-exposures state, the auto- 



covariance as given by Keane and Adrian (1990) is 



^ i JS)) = ^c(S) + ^ f (S) + ^p(5) + ^ d -(5) + ^ d+ (5). (3.5.41) 



Here the displacement-correlation peaks appears as 3^ D - or ^, D + (Figure 
(3.ii|). This is not a surprising result as the single-frame does not store any 



temporal information about the super-positioned realisations which causes a di- 
rectional ambiguity. 

There are techniques which can be employed to eliminate the directional 



ambiguity. A common method introduced by Adrian (1986^) is referred to as 
image-shifting, where the first realisation is spatially pre-shifted. The pre-shifting 
is performed by physically displacing the camera by a certain value during the 
two realisations. This solution also removes the pedestal which may make the 
extraction of some small displacements hard or impossible. 
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-^ 




Figure 3.11: Typical correlation plane for single-frame double-exposure mode. 



3.5.7 Estimation Values 

In the spatial correlation theory presented so far, it has been assume that the 



statistical properties of the PIV data are known. In particular, Equation (3.5.34) 
expresses a cross-covariance value which is calculated by some presumably known 
expectation values: 

^•uU^-i) = cov{I l (X ),!,(?+?£„)} 



= <E{ (l,(?) - £{I,(*)}) (lt(?+ ?? ;)l ) - E{I t (T+ T ?;)l )}) }. 

(3-5-42.) 
In this equation, the expectation operators perform ensemble averaging; that is, 
averaging over all possible values of the argument. In practice, when observing 
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a process, there is no pre-knowledge of that process available to describe the fea- 
tures of the process; however, it is possible to 'estimate' the features of the process 
by the observation itself. This is the fundamentals of the statistical inference tech- 
nique ( Hogg and Tanis|i996 ). 



Let the mean of the observed values be called the observation averaging. It 
can be shown ( Priestly 1981) ) the estimated values of the features of a static 



ergodic H process can be obtained when the ensemble averaging is substituted 
by the observation averaging. Implementing this to PIV data, the associated 
cross-covariance can be estimated when the ensemble averaging is replaced by 
the spatial averaging. The final form of the estimated cross-covariance is given 
later in § [3.6.2.1 , 



3.6 Digital Analysis 

Traditionally, the spatial correlation analysis method was applied to PIV data 
using hardware-based techniques. These included semi-optical methods, in- 



troduced initially for LSV data analysis by Burch and Tokarski (1968) using 
Young's fringes, and fully-optical methods, as suggested by Coupland and Hal 
liwell (E988J). The disadvantages of the opto-mechanical methods, like com 
plexity and temporal consumption, caused it to become obsolete and to be sub- 
stituted by software-based methods which have been increasingly efficient and 
economical. 

This section contains an introduction to the digital interface to the spatial 
correlation theory which, in a more general form, is known as Digital Particle- 
Image Velocimetry (DPIV). 

3.6.1 Pre- Processing 

Before the correlation-based digital image evaluation, some verification of the 
PIV data suitability may be necessary. This includes checking whether there is 
enough data to produce valid output and mapping the data to the digital space. 
The PIV data may be directly digitally recorded, or it may be digitised after 
recording. In both cases the digitisation parameters must be chosen in such 



I4 A process is said to be ergodic when, for infinite number of observations, the observation 
averaging is the same as the ensemble averaging. 
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a way that the data loss is optimised. The following subsections discuss these 
aspects of PIV data pre-processing in more detail. 

3. 6.1. i Image Sampling 

In the first step of digitising the particle-image field, the data space is divided into 
discrete components named pixels. With the sampling impulse response defined 
as 

-> f 1/aL if |x| < A px /2 A |£| < A px /2 

:h"i a iU)= p , (3.6.1) 

I elsewhere 

with A px being the size of the pixels, the discrete irradiance field can be obtained 
from the continuous irradiance field: 

Ia(Xa) = {JCi a i © I}(A px ."x A ), (3-6.2) 

where x" A = (x A ,y A ) is a dimensionless discrete vector. Using this definition, the 
associated auto- and cross-covariances can be derived as 

^i a (^a) ={(^i a i©:Hi a i)©^i} [a px .(7 a )] , (3.6.3) 

and 

^i A i A (^ A ) = h(z.7 n ){m l&l ®j{ lAl )®^ l }^A px .{7 A )- i^ , (3.6.4) 



(Westerweel 1993 §3.2), with ~? A being the dimensionless discrete correlation 



space. 

3.6.1.2 Data Bandwidth 

The bandwidth of a two-dimensional matrix with elements M.^ is defined as the 
value of |i — )\ such that M t j 7^ 0. In the case where the bandwidth of a matrix 
is a finite value, that matrix is said to be band-limited. According to the sampling 
theorem, a band-limited matrix is reconstructible if it is sampled at a minimum 
rate of twice of its bandwidth, also known as the Nyquist rate. 

Assume that the Fourier transform of the irradiance field is given by the 
matrix $ [I] where # denotes the Fourier transform operator. It can be shown 
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that # [I] is 'nearly' band-limited (Priestly 1981I. By sampling PIV data, the 



irradiance field may be treated as a signal which is reconstructed in the discrete 
form. The sampling rate, at which the data can be recovered, may be determined 
by the bandwidth of # [I] . 



Westerweel (1993 , §3.5) gives the bandwidth of the PIV imaging system as 



W P 



0.86 



(3.6.5) 



where fA is the Airy disk radius given by Equation (3.4.24I. Since the radius of 
the particle-image, as expressed by Equation ( 3. 4.27ft is larger than f A , Wp acts 
like an upper limit for the bandwidth of J [I], which is the PIV data bandwidth. 
According to this, for a typical imaging system, a data sampling rate of 64 mm -1 
is quite sufficient. 



3.6.1.3 Quantisation 

While data sampling converts the argument of the irradiance l[x) to a discrete 
value, quantisation digitises the value of the irradiance itself. The digitisation is 
a mapping from real values to some limited integers, also known as number of 
colors. 

Quantisation is an irreversible mapping which adds some noise to the output 
data. It has been shown ( W esterweel|i993 , §3.3) that the noise produced by an 
8-bit quantiser is negligible. That is, digitising the irradiance field to an image of 
256 colors — which is not necessarily, but typically, in gray-scale — preserves the 
data for processing. 



3.6.1.4 Valid-Data Yield 

An inappropriate value of the effective particle-image density C 01 > described 
by Equation ( 3.5.39) ), can cause a noise peak in the correlation plane which is 
stronger than the displacement peak. Introducing the quantity valid-data yield 
helps to optimise the value of C 01 • 



Similar to the discussion in § 3.5.1 , the probability of finding exactly k particles 



in an interrogation area may be written as 



N 



^k; 



X, 



OI e -X i 



k! 



(3.6.6) 
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where IN oi> which could be non-integer, is the effective number of particle-image 
pairs per interrogation area given by 



IN 



C 01 Si 



01 — c 01 j I 



(3-6.7) 



Assume k' is the number of particle-image pairs per interrogation area that is 
sufficient to remove the majority of the noise peak in the correlation plane. The 
valid-data yield is defined as the probability of having k' or more particle-image 
pairs in an interrogation area and is given by IJ 



k' 



IV (IN 01) =P k > k , 



:n. 



1 ~Lpk 



■N, 



(3.6.8) 



k=0 



Figure (3.12) shows some different values of F k /(INoi). It can be seen when the 



o.s 



n.o 



0.4 




Figure 3.12: Some values of valid-data yield. 



effective number of particle-image pairs per interrogation area is approximately 
between 10 and 15, the probability of finding at least 4 particle-image pairs per 
interrogation area (k' = 4) is almost 1. Choosing IN 01 to be a larger number 
causes overlapping of particle-images and should be avoided. 



I5 Westerweel I1993I § 3.6.4) gives the valid-data yield, in our notation, as T^'(Co\). This 
cannot be correct as the argument of the valid-data yield function must be a dimensionless num- 
ber. 
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To summarise, setting the interrogation area to have, roughly, 10-20 particle- 
image pairs is necessary to reduce the noise in the correlation plane. 



3.6.2 Data Processing 
3.6.2.1 Implementing Estimation 



Given the discrete irradiance field described by Equation (3.6.2) and following 



the arguments presented in §3.5.7 in this section the estimation value of the 
cross-covariance of discrete PIV data is presented. 

Consider two interrogation areas I a (Xa) and I a '(Xa) with the same size of 
N x N pixels. Given that the spatial average of discrete irradiance is 



N N 



i A = (i A (x" A )) = ^ Y Y MjTa). 



(3-6.9) 



it can be seen that the the mean value is the same for all realisations: 

<li(x A )> = (li , (x A )) = I A . 



(3.6.10) 



The cross-covariance of the two realizations I A (x" A ) and I a (Xa) over a shift 
~?" A can be estimated as 



N N 



£i A i A (7 A ) = ^2 X. Y- [^(^aI-IaJ [Ia(x a + 7 a ) -i a J , (3.6.11) 



x A = l y A = l 



where subtracting the mean values removes the background noise values %, c and 
% ¥ discussed in § 3.5.6 (We sterweel 1997} . The estimated cross-correlation is 
given by 



{I A *I A '}(7 A ; 



^-I A I^( S a] 



where the sign * denotes the correlation operator. 



(3.6.12) 



3.6.2.2 DFT-based Correlation 

The digital calculation of the measured discrete correlation plane, given by Equa- 



tion (3. 6. 11), is computationally intensive and, in practice, when the calculation 
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is needed to be repeated many times, is not functional. One way to overcome 
this problem is to perform the computations in the frequency domain, where the 
processing is less expensive, by making use of Discrete Fourier Transforms {DFT). 



Using the cross-correlation theorem, Equation (3.6.12 1 can be alternatively writ- 
ten as 

{ii*0(7 A ) = %i- 1 [3#r[ii] s^rai']] (7 A ), (3.6.13) 

where T>$Z denotes the discrete fourier transform. 

Usually, the DFT algorithm used in computations is the classical Fast Fourier 
Transform (FFT) method introduced by Cooley and Tukey (1965 1. The al- 



gorithm reduces the processing for an interrogation area of N x N pixels by a 
factor of 21og 2 (N)/N 2 . One disadvantage of the classical FFT is that it is only 
efficient when N is a power of 2, imposing restrictions on the interrogation area 



size. However, there are many other alternative algorithms (Press et al 



1992 



Chapter 12) which can lessen the restrictions and increase the interrogation area 
size choice. 



3.6.2.3 Sub-Pixel Peak Detection 

The information stored in the discrete correlation plane, given by {I A * I A }(T A ), 
allows the position of the displacement peak to be determined with a resolution 
of less than a pixel. In fact, since the correlation peak always covers more than 
one discrete unit in the correlation plane ( Westerweel|i997 ), the displacement 
peak can be detected with sub-pixel accuracy. Based on this idea, Willert and 
Gharib ( 1991} introduced the sub-pixel peak detection theory by interpolation 
which is essential for calculating second order quantities, e.g. vorticity or pres 
sure. 
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Chapter ^j- 

Acoustic Velocimetry 
Using Particle-Imaging 



Everything you can 
imagine is real. 

Pablo Picasso 

Laser velocimetry methods are ordinarily applied to fluid flows under no 
or negligible external forces. Consequently, when implemented on sound 
fields, PIV principles should be validated in order to meet the new condi- 
tions. This chapter examines aspects of PIV theory and discusses the limitations 
of PIV when applied to sound fields. 

4.1 Acoustic Velocimetry Techniques 

Unlike the scalar quantities associated with fluids in sound fields, particularly 
pressure and temperature, measurement of acoustic velocity vectors is not straight- 



forward. Figure (4.1 1 shows the most common acoustic velocimetry techniques 
classified in terms of their level of intrusion during a measurement. The tech- 
niques are described briefly in this section. 

4. 1. 1 Intrusive Techniques 

In the late 19 th century, Rayleigh described a system for measuring acoustic ve- 



locity (Rayleigh 1945, §253b). After observing the torque that a non-symmetrical 
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Acoustic Velocimetry 



LDV 
(Taylor ipy6) 



r 




Rayleigh Disc 
(Rayleigh ip4s) 



Two-Microphone 
(Fahy ipyy) 



HWA 
(Pratt el al. ipyy) 



Microflown 
(van der Eerden et al. ipp8) 




FDV 
(Buick et al. 2002) 



Figure 4.1: Acoustic velocimetry techniques are grouped according to intrusiveness. 



obstacle of appropriate dimensions experiences in a sound field, Rayleigh sug- 
gested using the deviation of a disc floating in a fluid from its equilibrium po- 
sition to calculate the acoustic velocity of the fluid. The technique, which is 
referred to as the Rayleigh disc, has now been consigned to history and is no 
longer in use. 

One way to measure acoustic velocity vectors is to calculate them from scalar 
quantities which can be measured more easily. For example, a two-microphone 
method may be applied to extract velocity from pressure ( Fahy|i977 ). 



Thermo-electric methods like Hot-Wire Anemometry {HWA) (Pratt et al. 



1977) or Hot-Film Anemometry (HFA) offer a better accuracy at the expense of a 
more complicated set-up. The Microflown (van per Eerden et ALJ1998 1, which 
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is based on a thermo-electric structure, is convenient to use but has an upper 
limit for velocity measurement which is usually lower than the range of interest 
in acoustic studies. 

Because of their physical presence in a sound field, and because their dimen- 
sions are not negligible when compared to the wavelength of the sound field, all 
of the devices used in the techniques mentioned so far have the disadvantage of 
being intrusive in nature. 

4.1.2 Non-Intrusive Techniques 

As they are optically based, the laser velocimetry techniques inherit the charac- 
teristic of non-, or at least low-, intrusion. Of these, LDV (LDA) and PLV are of 
special interest. 

Taylor ( 1976} was the first to apply LDV to resonant tubes. The technique 



involves measuring acoustic velocity by the Doppler shift of waves scattered 
from tracer particles within small volume found by the intersection of two laser 
beams. The technique is now well-established and widely used; however, it suf- 
fers from the lack of a whole-field feature. 

The use of PLV in the study of sound fields has only recently taken off due to 
previous limitations in capturing images of high-speed flows. The applications 
of this technique in sound fields may be grouped according to the type of fluid 
markers used in the measurements. 

Fluorescent Dye Velocimetry (FDV) is a PLV-based technique in which dye is 
used as the molecular marking (§3.2.1) for fluid visualisation (Buicket al.|2qo2J). 



The method has been applied successfully to acoustic streaming measurements. 
In comparison to PIV, FDV has the advantage of eliminating the interaction 



between sound field and particulate fluid markers (§4.4.2); however, there are 
limitations in applying dye to the fluid under investigation. 

Using particulate marking in PLV sound field measurements actually involves 
acoustic PIV, which is the focus of this chapter. 

4.2 Acoustic PIV Applicability 

Investigation of the limitations of applying PIV to sound fields may be imple- 
mented in two parts. First, the restrictions of PIV as an existing technique is 
presented. This includes positioning the spatial and temporal ranges of PIV in 
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the framework of sound fields. Later, in §4.4 the effects of the presence of sound 
fields in PIV measurements are discussed. 



4.2.1 PIV Ranges 



A detailed description of spatial resolution of particle-images is given in §3.4 At 
the time of this study, a typical CCD camera has a resolution of the order of 
1000 x 1000 pixels together with the pixel spacing A R < 10 um. Given some 
rough values for the laser wavelength A ~ 500 nm, magnification M « 2, f- 



number f # ss 1 and applying Equation (3.4.27 1 it may be seen that a 10 urn 
displacement is mapped to almost 3 pixels. The associated uncertainty is ap- 
proximately 1 um. 

For improving the spatial resolution of PIV down to micron and sub-micron 
level, some substantial work has been carried out leading to a method known 



as Micron-resolution PIV {micro-PIV or \xPIV) introduced by Meinhart et al. 



(1999). Their uPIV set-up could map a 200 nm displacement to almost 3 pixels 
(spaced by 6.8 um) with an uncertainty of about 53 nm. The technique requires 
sub-micron seeding particles for tracing the flow faithfully and sub-micron width 
of the illuminating light sheet for avoiding out-of-the-plane motions. 

The temporal resolution of a PIV set-up depends, at least, on the following 
parameters: 

Illumination. Intensity of illumination, duration of each illumination 
and time separation between illuminations. 

Recording. Sensitivity of recording medium, duration of each recording 
and time separation between recordings. 

Synchronisation. The duration of each illumination and each recording 
together with the corresponding time separations must be matched. 

Of these, the right temporal parameters of the recording parameters seems to 
be harder to achieve. Recent advances in digital imaging has eliminated a big 
part of this problem by the use of interline CCD chips, allowing recordings with 
a time separation of the order of nano-seconds. However, the technology still 
suffers from a lack of flexibility in recording durations. This deficiency may 
be covered by taking advantage of the more flexibility available in the duration 

56 



4-2. Acoustic PIV Applicability 



of illuminations. The limitation of this method is determined by the trade- 
off between the flexibility of illumination duration and illumination intensity; 
that is, shorter illumination durations may be reached at the expense of lower 
illumination intensity. 

4.2.2 Sound Field Ranges 

One feature of acoustical phenomena is that their characteristic quantities have a 
broad range. In particular, this may involve the orders of 1 Hz to 10 8 Hz in the 
frequency domain while, in different media and depending on the pressure level, 
the acoustic displacement of interest could vary from 10 -8 m to 10 2 m. This 
makes the development of a versatile measurement system a challenging task. 

Sound fields in the linear regime simply involve acoustic motion, that is, the 
motion of the fluid in the field is directly generated by the sonic source. As the 
Reynolds number of the fluid — identified by the ratio of the inertial force to 
the viscous force — gets higher values, the fluid transits to the non-linear regime. 
The fluid in the non-linear regime may be purely in the acoustic motion phase, 
as when energy is transferred to the higher modes, or it could transit to the hy- 
drodynamic motion where secondary flows are formed as a result of time-averaged 
fluctuations with a non-zero mean. These hydrodynamic motions are mainly 
known as acoustic streaming (NYBORGI1965; Lighthill|i978}) or steady streaming. 1 



Figure (4.2 1 illustrates sound field ranges in the frequency and space domains 
in comparison with the PIV ranges discussed previously. The sound field ranges 
are grouped into acoustic and hydrodynamic motions and the ranges of each 
group, as explained later, are presented based on the corresponding characteristic 
parameters. It is assumed that the PIV measurements sample the sound field 
under investigation at a rate of ~ 10 in space and time. Hence, as shown in 
the figure, the 'applicable' PIV ranges are smaller than their actual values by one 
order. 

The characteristic frequency and displacement of acoustic motions is essen- 
tially the frequency and amplitude of the associated oscillatory motions. In 
air, typical amplitudes of linear oscillatory motions with 100 to 200 dB SPL 
((P re f = 20 uPa) and a frequency range of 1 to 2 kHz is reported to be of the 



'Riley (2001I suggests the term 'steady streaming' as an alternative to 'acoustic streaming'. 



This is due to the fact that the associated phenomenon involves incompressible flows, while, 
'acoustic streaming' implements compressibility of the fluid. 
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Figure 4.2: A rough illustration of acoustic PIV ranges based on the characteristic 
frequency and displacement of the fluid. The product of the characteristic parameters 
at each point gives the corresponding velocity amplitude. It is assumed that the PIV 
(\iPIV) measurements sample the fluid motion at a rate of ~~ 10 and the cut-off 
velocity of the fluid is around 10 mis. The dashed lines represent the PIV (\iPIV) 
borders. 



order of sub-millimeter ( Hann and Greated i997afb ). Thus, acoustic motion 



is, at least in some regions, measurable by PIV. In contrast, based on Vignola 



et al. 



of 



(1991), amplitudes in water with 180 dB SPL (T ref = 1 uPa) at a frequency 
2 kHz are of much lower sub-micron values which are quite out of the 
PIV borders. Nevertheless, it seems to be possible potentially to apply uPIV to 
acoustic motions in liquids. It should be noted that applying uPIV to sound 
flow passages of dimensions much greater than the micron order cannot yield 
relatively whole-field results, being closer to a point measurement technique. 

Increasing the sound intensity level in order to gain higher amplitudes of 
acoustic motions may cause two undesirable effects. First, this increments the 
acoustic velocity which demands higher temporal PIV sampling and, conse- 
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quently, shrinks the PIV ranges. Secondly, increasing the amplitude in this way 
can create some unwanted non-linear events such as distortion or streaming. 

In hydrodynamic motions, the characteristic parameters are a few orders 
greater than the ones of acoustic motions. Acoustic streaming velocities are 
normally of the order of < 1 cm/s and rarely up to a few meters per second 
(Akulichev 1971). This well-fits the hydrodynamic motions within the PIV 



ranges. 



4.3 Previous Studies On Acoustic PIV 

Based on the discussed acoustic ranges, the applications of PIV to sound fields 
can be mainly divided into the measurement of acoustic and hydrodynamic mo- 
tions. The following is an overview of these applications where some of them are 



also covered in Campbell et al. (2000) 



4.3.1 Acoustic Motions 

Due to the restrictions in applying PIV to acoustic motions, namely difficulties 
associated with capturing images rapidly enough to track fast tracer particles, the 
literature on this topic is quite sparse. 

In the first attempt to use PIV to study a sound field, due to the lack of 
illumination intensity for short exposures and/or the difficulty in sending illu- 
mination pulses fast enough, the image recording process was carried out over 
a few acoustic cycles (Hann 19951 Hann and Greated i997a[b] >. The acquired 
data was a set of streaks rather than particle-images and consequently required a 
slightly different data analysis method to that used in general PIV. As the image 
recoding process takes more than half of an acoustic cycle, no phase information 
can be obtained using this method. 

The direct measurement of acoustic displacement in air using real particle- 



imaging was introduced by Blackshire ( 1997 1 and later extended by Humphreys Jr 



et AL. 



( 1998] ). The output of the measurement is the 'absolute displacement' of 
fluid particles from their equilibrium which is not the same as the 'instanta- 
neous displacement' that leads to velocimetry The technique is based on two 
independent laser sources, with a static delay of half an acoustic cycle, which 
are synchronised with the sound field by a microphone trigger. The displace- 
ment data is extracted from a single-frame double-exposure recording followed 
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by auto-correlation analysis. The advantage of this system is that it extracts phase 
information as well as the amplitude and direction of acoustic displacement. 



4.3.2 Hydrodynamic Motions 

The application of PIV to the study of different types of acoustic streaming, that 
involve velocities which are much lower than the velocities associated with typical 
acoustic oscillatory motions, demands less challenging conditions. 

Initiated by Sharpe et al. (1988k a considerable amount of research has 



been undertaken on PIV applied to hydrodynamic motions; namely, Rayleigh 
acoustic streaming in air (Hann and Greated 1997c! Rockliff et al 



2000 



Rockliff and Greated 2ooi[ Skulina et al. 2003} Joflansson et al^|2oo3| ) 



and Eckart acoustic streaming in water ( Cosgrove et al.|2qoi >. Many of these 
studies involve double-frame single-exposure recordings with an exposure time 
of the order of a few tens of milliseconds or more. 



4.3.3 Novelty 

The use of PIV in the study of sound fields has provided a better qualitative 
understanding of acoustic behaviour, particularly, in complex geometries. In 
addition, quantitative results such as standing wave pattern and the impedance 
of low-speed flow passages have been achieved. 

However, from the flow visualisation point of view, not much remarkable 



progress can be seen in recent works. Figure (4.3) compares recent visualisations 
with some similar cases performed in the 1930s. The lack of progress can be 
seen particularly in the measurement of high-speed oscillatory motion where no 
substantial work has been reported on the capturing of particle-images within 
one acoustic cycle. 



4.4 Effect of Sound Fields On PIV Measurements 

In its general form, PIV is known as a non-intrusive technique, as has been 
well-verified through applications such as the measurement of turbulent flows. 
However, the presence of a sound field in the fluid under investigation introduces 
various new phenomena which may alter the fundamental assumptions that PIV 
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a) Brandt and Hiedemann (1936). 






c) Andrade \1931) 



d) Rockliff (2002) 



Figure 4.3: A comparison between flow visualisation on PIV measurements in sound 
fields in the ipjos to similar studies in recent research. Figures (a) and (b) show 
the visualisation of acoustic oscillatory motion over a few cycles. Figures (c) and (d) 
correspond to visualisation of hydro dynamic motion generated by Rayleigh streaming. 



provides a non-intrusive measurement (§3.21. Accordingly, a careful considera 



tion of any acoustical effect on the laser velocimetry configuration must be taken 
into account. 

To study the level of intrusion incurred during acoustic PIV measurements, 
it is convenient to decompose a PIV measurement into sound field, light field 



and tracer particles components, as illustrated in Figure (4.4 1. Under this as 



sumption, the measurement can be considered non-intrusive provided that none 
of these components interact with each other or if the interaction is negligible. 
The amplitudes of the interactions depend crucially on the characteristics of the 
sound field, like frequency and sound intensity, and the material and optical fea- 
tures of the fluid and tracer particles. Hence, the amount of intrusion during an 
acoustic PIV measurement is very specific to that measurement. The following is 
a discussion of the common interactions which stop acoustic PIV measurements 
from being completely non-intrusive or providing reliable output data. 
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-I Tracer Particles J- - ' 




Figure 4.4: The intrusion parameters in acoustic velocimetry using particle-imaging. 



4.4.1 Sound Field - Light Field Interactions 

The presence of sound waves in a medium is a source of optical inhomogeneity 
which makes the light scattered from an object in the medium experience a gra- 
dient of the refractive index. This causes the apparent optical properties of the 
object, such as the position, to differ from the exact values. 

In the case of ultra-sonic waves, the interaction between sound and light fields 



can be explained by acousto-optics (Sapriel 1979). In acoustic LDV measure 



ments, the acousto-optical effect has been studied for propagating (Jack et al. 



1998 1 and standing sound waves (Jack et al. 1999) in water. It has been shown 
that sound waves may intrude the Doppler signal in a non-negligible way. 

The aero-optical effect on the distortion of PIV data has been pointed out 



by Elsinga et al. (2004 1, reporting three types of errors introduced by the op- 



tical inhomogeneity of the medium under investigation: the position and veloc- 
ity errors and particle-image blur. The first two errors are systematic while the 
last one increases the signal-to-noise ratio of the PIV data. In addition, it was 
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demonstrated that the systematic errors are related to the second derivative of the 
refractive index in the velocity direction. 

In acoustic PIV, the instantaneous refractive index of the fluid under obser- 
vation is proportional to the acoustic pressure. According to this and based on 



the results by Elsinga et al. (2004 1, the source of the velocity error is only sig- 
nificant when the second order pressure gradient in the direction of velocity is 
remarkable. 

4.4.2 Sound Field - Tracer Particles Interactions 

The interaction between sound field and tracer particles in a PIV measurement 
is mutual and may have a strong intrusive participation. 

The presence of tracer particles in a fluid can cause attenuation of the prop- 



agation of sound in that fluid (Kinsler et al. 2000 §8.10). In fact, the more 



faithfully the tracers follow the fluid, implying that more momentum from the 
sound waves is transferred to them, the greater the attenuation of the sound 
wave. This leads to a sort of 'uncertainty' between the accuracy of the acoustic 
PIV measurement and the disturbance of the measured acoustic wave. 

Conversely, the effect of sound fields on tracer particles, beyond helping to 
measure the fluid motion by forcing the tracers to follow it, takes the unwanted 
form of spatial redistribution of the tracer particles. This happens as a result of 
migration of the particles introduced by a non-zero averaging motion. Moreover, 
the different motions of the tracer particles may cause them to collide and adhere 
to each other which leads to a size redistribution. As will be shown later, a con- 
sequence of the size redistribution is the velocity redistribution of the particles. 

The remainder of this chapter focusses upon the sound field - tracer particles 
interactions. 

4.5 Motion of Small Particles in Standing Sound 
Fields 

4.5.1 Terminology 

The term aerosol refers to a suspension of fine solid or liquid particles in a gaseous 
medium. A similar suspension in a liquid is referred to as a hydrosol. The scope 
of this study is limited to aerosols. 
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In the presence of sound fields, rigid aerosol particles are defined as volumes 
with no sound waves inside them, while sound waves can exist inside deformable 
aerosol particles. 

The motion of fluid around an aerosol particle applies a force to it. The 
elements of this force which are parallel and perpendicular to the fluid motion 
are known as the drag and lift forces, respectively. When the fluid motion is 
oscillatory, as in an acoustic field, the temporal average force that aerosol particles 
experience is named the drift force. 

The Reynolds number of the aerosol particles helps to describe how the fluid 
flows around the particles. The number is defined as the ratio of the inertial and 
the viscous forces which, for a spherical particle, may be given as 

2r D I v"— v" D | 

:N Re = ^ p -, (4.5-1) 

V 

where r p is the radius of the particle, | v"— V p | is the relative velocity of the particle 
and the fluid and the kinematic viscosity coefficient is denoted by v. 



4.5.2 Equation of Motion 

In a turbulent flow, the motion of spherical particles with dimensions much 
smaller than the minimum wavelength of the turbulence may be described by 



the Basset-Boussinesq-Oseen (BBO) equation (Boussinesq 1885 ; Basset|i888 



Oseen|i927L A general form of the BBO integro-differential equation, which is 



also valid for the spacial case of laminar flow, can be written as 

An 3 d _* An 3 ->■ 

T iPp S v p = - T r3vp 

1 An , d ,_+ .y , 

+ 2T r p po dl (v - v f ) 

An a 1 ,_+ _,. . 

+ T r p pp ^ (v " vp) 



+ 6r; v /7tp ri 



v - vJdt' 



t0 x/t-t'dt' 



v> 



+ Fc» (4-5- 2 ) 



64 



4-5- Motion of Small Particles in Standing Sound Fields 



where p p is the density of the aerosol particle and r\ is the dynamic viscosity of 
the fluid. 

The first term on the right hand side of this equation describes the force 
generated by the pressure gradient Vp in the fluid around the particle. 

The second term is the force due to the apparent mass of the particle relative 
to the fluid. This force is the outcome of the direct exchange of momentum 



between the particle and the fluid, assuming that the fluid is ideal ( Lamb||i932 ). 

The third term of the equation is due to the viscous drag force with 1/t d as 
a time constant for the momentum transfer given by 



1 3 po 

t d 8 P P r i 



c d:t-Hv - v p |, (4.5.3) 



where Cd is the drag coefficient. 

The fourth term of the BBO equation (4.5.2) is the Basset history integral 
force. This force is due to the non-uniform relative motion of the viscous fluid 
depending crucially on the path of the motion of the particle. 

The last term is due to any potential field force such as gravity or electromag- 
netic forces. 



4.5.2.1 Assumptions 

The BBO equation is valid for a wide range of parameters. The scope of interest 
throughout this study is defined by the following assumptions. 

Assumption 4a. The aerosol particles have a spherical shape and they are rigid 
in sound waves. 

Assumption 4.H. The particles do not have any kind of interaction with each 
other. 

Assumption 4.iii. The dimension of the aerosol particles is much smaller than 
the wavelength of the sound field, 

?«1- (4.54) 
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Assumption 4.1V. The density of the aerosol particles is much greater than the 
equilibrium density of the host fluid, 



Pp 
Po 



»1. 



(4.5-5) 



Assumption 4.V. The velocity of the aerosol particles relative to the host fluid 
changes slowly so that 

d ,- - 



dT lV - V r J 



« 1. 



(4-5.6) 



Assumption 4.V1. The aerosol particles have low Reynolds numbers implying 
that the drag coefficient is given by the Stokes' law regime as (Stokes|i847> 



C 



24 



D 



DV 



Re 



(4.5.7) 



4.5.2.2 Equation of Motion in Standing Sound Fields 



The preceding assumptions enable the BBO equation (4.5.2 1 to be greatly simpli- 
fied. For example, taking into account Assumption ( 4.iv| , the pressure gradient 
and apparent mass terms are negligible. In addition, Assumption (4.V) ensures 



that the Basset history integral value can be ignored. Finally, Assumption <4.vi) 



determines the final form of the viscous drag force giving the one-dimensional 
equation of motion of aerosol particles as 



An , d 
yr p P P 



dt^ 



6nr\ r p (v - v p 



(4.5.8) 



Introducing the relaxation time of an aerosol particle as a function of the size of 
the particle: 



TplTpj 



2 Pp Tj 

9 po y' 
the equation of motion can be re-arranged to 



(4-5-9) 



'dt 



(4-5.10) 
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In the particular case of a standing wave sound field, substituting Equation 



(2.4.9 1 into the equation of motion yields 



d. 
'dt' 



T p — v p + v p = Vq sin(kx) cos(cut) 



(4-5-n) 



The general form of this differential equation may be given as (Polyanin and 
Zaitsev||i995| 



v p (x, t) = v x Up sin(cut- cp) +v x \i & e 



-t/T„ 



(4.5.I2a) 



v x = Vq sin(kx) 



(4.5.12b) 



<P 



tan 



, tDT n 



(4.5.12c) 



Up = cos(cp) 



(4.5.12c!) 



u g = sin((pj, 



(4.5.12c) 



where cp is a phase shift angle and u p and \i % are the entrainment coefficient and 
the flow-around coefficient, respectively. After a time approximately equal to T p , 
the second term in Equation ( 4.5.12a! is negligible and the particle is approx- 
imately in a steady-state. Considering this, the velocity of an aerosol particle 
in a sinusoidal standing wave is sinusoidal in nature with the spatial amplitude 
reduced by a factor of cos ((p) and a temporal delay of cp introduced. 



4.5.3 Velocity Distribution 

In practice, an ensemble of aerosol particles has a size distribution. According to 
Equation ( 4.5.12a! , the velocity of an aerosol particle in a sound field is a function 
of its size. This implies that, at any space-time, aerosol particles have a distribu- 
tion in the velocity space. One immediate consequence of this argument is that 
the number of tracer particles which follow the fluid faithfully is a function of 
space-time. Consequently, when performing an acoustic PIV measurement, a 
pre-existing knowledge of the velocity distribution behaviour of the tracer parti- 
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Figure 4.$: Values of the entrainment coefficient for different frequencies. 



cles in the sound field under investigation may be crucial. 

4.5.3.1 Fluctuation of Number of Illuminated Particles 

The tracer particles which contribute to velocimetry, distinguished by illumina- 
tion, are a sub-set of the whole ensemble of particles. Given a fixed illumination 
volume, in the presence of a sound field, the number of contributing particles is 
a function of the dimensions and position of the visualisation volume, as well as 
of time. The aim is to calculate this number at any space-time as a pre-requisite 
for obtaining the velocity distribution. 

Consider a resonant tube with aerosol particles homogenously distributed in- 
side it. It is assumed that the velocimetry measurement is performed in a narrow- 
enough band, with a length of Axl, across the tube so that the velocity within 
the illumination volume has an almost constant value implying the criterion 



Ax L < A. 



(4-5-I3) 



Let the instantaneous number density of the aerosol particles be n p ; that is, the 
number of particles per unit volume at any time. The aerosol particles and the 
host fluid can be treated as a multi-phase system where each phase fulfills the 
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continuity condition independently. Hence, the equation of continuity for the 
aerosol particles may be written as 



3 t Tip + 3 x (n p Vp) = 0, 



(4-5-I4) 



Expressing the instantaneous number density as equilibrium number density n Po 
and fluctuating number density n p : 



n p = n p + 6n p> 



the equation of continuity can be linearised to 

d t &n p + Up o 9 x v p = 0. 



(4-5-I5) 



(4.5.16) 



Using the boundary condition of Equations (2.4.4) anci (2.4.6), a solution of the 
linear continuity equation is 



6u p = — n Po — cos(cp) cos(kx) cos(a>t— cp), 



(4-5-I7) 



giving the instantaneous number density as 



Tip(x,t) 



n Po ( 1 — n Po — cos(cp) cos(kx) cos(cut— cp) J . (4-5-i8) 



The number of illuminated particles in the volume (Azl WAxl) can be given by 

fX+Ax L 



^ p = Az L W 



TLp(x) dx, 



(4-5-I9) 



where W and Azl are the width of the tube and the thickness of the light sheet, 
respectively. After integration, the number of illuminated particles may be ex- 
pressed in terms of equilibrium and fluctuating parts by 



w. 



here 



^p = ^p o + 6^ p , 



^ Po = (Az L Ax L W)n, 



Po' 



(4.5.20a) 



(4.5.20b) 
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and 

Kn , r ~ r v , , . . . /sin(kx-kAx L )-sin(kx) , 

63V p = 3V Po — cos((p) cos[cut- cp) I — ] . (4.5.20c) 



Using Equation (4.5.13}, when kAxL <C 1, the fluctuation in the number of 



illuminated particles can be approximated as 



6ZNT p = N Po — cos(cp) cos(cut — cp) cos(kx). (4.5.21) 



Equation (4.5.21 1 shows that the value of 6N p depends mainly upon the ratio 
of the acoustic velocity to the phase velocity. In linear sound fields this ratio is a 
small number implying that the fluctuation is negligible in the linear regime and 

:n p = :n Po . (4-5-22) 

This value of the instantaneous number of illuminated particles will be used to 
calculate their velocity distribution. 



4.5.3.2 Calculation of Velocity Distribution 

As a consequence of homogenous spatial distribution, the illuminated aerosol 
particles in an acoustic PIV measurement have almost the same size distribu- 
tion as the whole ensemble. The probability density function (PDF) of the size 
distribution is defined by 

d^ Po =/ r (r p )dr p , (4-5- 2 3a) 

where dN Po represents the number of illuminated particles whose radii are be- 
tween r p and r p + dr p . The normalisation condition gives the total number of 
illuminated particles at any space-time which, following on from the discussion 



in §4.5.3.1, may be approximated to its equilibrium value: 

00 
/ r (r p )dr p = W Po . (4.5.2.3b) 

Commonly, the size distribution of aerosol particles is described by a log-normal 
distribution. 
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In a similar way, the PDF of the velocity distribution of the illuminated par- 
ticles and the associated normalisation condition may be written as 



cW 



Pn 



:/ v (Vp)dV p 



(4.5.24a) 



r+ v o 



/ v (vp)dv p 



N 



Po 4 



(4.5.24b) 



-v 



Combining Equations (4.5.23a) and (4.5.24a) and expressing each differential el- 
ement of the velocity space in terms of a differential element in the size space, 
the PDFs of velocity and size distributions are related to each other by 



/v(v p ) a^Vpa^cpdrp =/ r (r p )dr p 



(4-5-2-5) 



Typical values for the phase shift cp of the tracer particles commonly used in 
acoustic PIV measurements are much smaller than unity. Considering this, after 



some algebraic manipulation, Equation (4.5.25) reads 

/v(Vp) = 



2cpvo sin(kx) sin(cut — 2cp' 



/t(t p ), cp < 1. 



(4.5.26) 



Since the area under the distribution curve is assumed to be constant, a rough 
estimation of the width of the distribution may be given by 



cr Vo Vq |sin(kx) sin(tut — 2cp)|, 



w. 



here 



2cpW 



Po 



'Vq 



r p/r(r p 



(4.5.27a) 



(4.5.27b) 



Equation (4.5.27a) states that at the extremum points of the acoustic velocity 



the velocity distributions of the aerosols are wider than the velocity distributions 
where the acoustic velocity is zero. The ratio of these velocity distribution widths 
is given by tan(2cp). 
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Chapter S 



Instrumentations 



Dear fr lend, 

all theory is gray; 

And green 

the golden tree of life . 

Johann Wolfgang von Goethe, Faust 



The experimental set-up, designed for acoustic PIV measurements, comprises a 
set of hardware which is controlled by, and whose output is analysed by, self- 



written software (Figure 5.1 1. The measurement cycle is initiated by the software 
sending a series of input signals to the hardware. The signals contain the re- 
quired information for both running the experiment and performing the mea- 
surement. The measurement cycle continues as the hardware output is received 
and analysed by the software. 

Typically, acoustic PIV measurements need to be run temporally as well as 
spatially, requiring a large number of inputs and the analysis of a high volume 
of output. This makes automating the measurement cycle a crucial factor in 
speeding up the process and cutting down on human error. To this end, the 
software layer has been designed to allow the number of hardware components 
to be kept to a minimum, enabling a large number of PIV measurements to 
be made in rapid succession via automation of the system. Some other benefits 
include making the system more flexible, time-effective and cost-friendly 

In the following, the hardware and software components are discussed in 
detail. 
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Dataflo 



Physical, electrical 
or optical connection 



Figure $.1: Experimental design. 



5.1 Hardware 



The hardware configuration is illustrated in Figures (5.1 1 and (| 5 . 2,)) . The input 



signals, explained in §5.2, are sent to a data acquisition board (DAQ) housed 
within a PC. The DAQ performs two main jobs; generating voltages to trigger 
the other devices and synchronising or delaying the triggered events. In this way, 
the DAQ is used to control the loudspeaker, to create the desired phenomena 
inside the tube; and the laser and digital camera, to perform the measurements 
on the tube. The measurement process ends with the recorded images from the 
camera being downloaded to the PC via a frame-grabber card and stored in a 
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digital medium for processing later. The technical specifications of the devices 



are given in Table (5.1 1 




Laser Optical Sound Field 
Set-up Set-up 



Figure $.2: Experimental set-up: hardware configuration. 



5.1.1 The Tube 

The transparent tube used in the study is made of Perspex 1 material. The ther- 
mal conductivity of Perspex is about 16% that of glass, making it possible to 



'Perspex is the UK trade name for polymethyl-methacrylate (PMMA). 
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De 



Specifications 



DAQ 

Loudspeaker 

Laser 

light sheet generator 

Camera 
Camera lens A 
Camera lens B 



IOtech DaqBoard-2000, 16 bit, 200 kHz 

Laser Quantum, 532 nm, 250 mW, C W solid-state 
Three cylindrical lenses (Figure 3.6c!, 



Focal lengthes: +19 mm, +200 mm and +60 mm 
PCO DoubleShutter Sensicam 
Yf = 12.5 - 75 mm 
Yf = 55 - 110 mm 



Tube Perspex 

Length 62.0 cm 

Square cross-section 3.5 cm x 3.5 cm 

Wall thickness 0.5 cm 



Table 5.1: Specifications of hardware devices. 



thermally insulate the volume inside the tube more efficiently than glass. More- 
over, allowing about 92% of visible light, Perspex is a better choice over glass in 
transmitting light energy. 

One end of the tube is terminated in an acoustically rigid manner while the 
other end is, via an adaptor, connected to the loudspeaker. The tube length is 
62.0 cm with a square cross-section of dimensions 3.5 cm x 3.5 cm and a wall 
thickness of 0.5 cm. 

The choice of a cuboidal tube shape is in order to avoid the flare problem 
which is common in cylindrical tubes due to the unwanted light reflected from 
their curved walls. The area in PIV images which is affected by flare suffers a 
severe lack of data. In principle, it is conceivable to overcome the flare problem 
when the light scattered from the tracer particles has a different wavelength to the 
light reflected from the walls. 2 In this way, the flare light may be easily filtered. 

5.1.2 Tracer Particles 

During the measurement, the tube under investigation is seeded using smoke 
for the tracer particles. The smoke is generated by burning incense of the joss- 



stick type. Figure (5.3) shows the size distribution of typical joss-stick smoke 



2 This is possible by dying the tracer particles using fluorescence when, after absorbing the 
illumination light, they can emit a light with a different wavelength. 
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particles measured by an Electric Low Pressure Impactor (ELPI) (Jetter et al. 



2002| . The peak of the size distribution is located between 2.6 and 4.0 urn. 3 As 
generated from oil-based material, the density of the smoke particles is expected 
to be of the order of 1000 to 2000 Kg/m 3 . 



Following the discussions in §4.5.2.2, at sound field frequencies up to a few 
kilo-Hertz, sub-micron aerosols whose density is much higher than of air have an 



entrainment coefficient near unity (Figure 4.5). This implies that, at low-enough 



frequencies, the smoke particles are expected to follow the air flow faithfully. 

While incense smoke is the most cost-effective way of visualising air flows, 
there are alternative methods which provide better control over physical charac- 
teristics of the tracer particles. One example is using aerosol generators which 
create micron-sized particles from almost arbitrary liquids. The tracer particles 
generated by this device have a sharp peak of size distribution. Moreover, in order 
to get particle-images of better quality, it is possible to increase the intensity of 
the light reflected from the tracer particles with the right choice of source liquid. 



5.1.3 Light Sheet 

The optical set-up used to generate the light sheet, for illuminating the tracer 



particles, consists of three cylindrical lenses, as shown in Figure (3.6c) on page 
25, with focal lengths given in Table (5.1). The generated light sheet is able to 



illuminate a region of around 1 cm high across the width of the tube, thickness 
of the light sheet has a Gaussian profile with a width of less than 1 mm. 



5.1.4 Laser 



The option exists to run the laser in either a CW mode (§3.3.1! or to trigger it 
by an external source to produce a beam with an intensity varying as an arbitrary 
function of time. 



'Commonly, in order to fulfil health and safety regulations, burning incense is preferred to 
cigarette smoke. Measurements by Jetter et al. 1 2002 points out that incense smoke particles 



are of sub-micron size and may be as unhealthy as cigarette smoke. 
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Figure 5.3: Size distribution of smoke particles generated by burning typical joss-sticks 
with apeak between 2.6 and 4.0 \wi. The data is adopted from Jetter etal. (2002) 



5.I.5 



Camera 



The digital camera device used in the set-up may be triggered, by an external 
source, to perform two exposures: a dynamic exposure, whose duration can be 
controlled by the trigger source, and a static exposure which lasts for a constant 
time of the order of a few milliseconds. The end of the first exposure is separated 
from the beginning of the second exposure by a delay of the order of nanosec- 
onds. 

In a PIV set-up, the recording process may be controlled by either illumina- 
tion (the laser) or by exposure (the camera). In the first method, depending on 
the purpose of the experiment, the camera shutter is left open to expose one or 
more laser illuminations. In the second method, the laser is run in the CW mode 
and the exposures are performed by the camera at the desired times and dura- 
tions. Due to the fact that the second exposure of the camera is not controllable, 
it is only practical to use the first method in this experimental set-up. 
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Figure $.4: Experimental set-up: hardware configuration. 



5.2 Software: Input Signals 

The following signals are used to control the hardware configuration: 

sound wave signal. The loudspeaker is driven by the corresponding sig- 
nal for the sound field of interest. Throughout, this signal is referred to as 
the sound wave signal. 

illumination signal. The laser is triggered by the illumination signal to 
produce 'quasi-pulses'; that is, pulses which are generated as a result of 
applying a square-function to the laser running in the trigger mode. Quasi- 
pulses provide the flexibility of being able to be triggered at any desired 
time. However, since there is no energy storage prior to sending the quasi- 
pulse, the output light may not be intense enough for some measurements. 

exposure signal. A square-function signal, similar to the illumination 
signal, referred to as the exposure signal, is sent to the digital camera device. 



Following the discussions in §5.1.5: the recording process is controlled by 
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the illumination signal. Consequently, in order to ensure that each PIV 
image is captured at the instant when the tube is being illuminated by the 
light sheet, the recording signal must be carefully synchronised with the 
illumination signal. 

The illumination and exposure signals are in the form of high-TTL. 4 De- 
pending on the aim of the measurement, various combinations of the input sig- 
nals may be used. The illumination signals used in this study are of identical 
duration and each is covered by one exposure signal. This configuration pro- 



duces multiple (double) frames - single exposure data (§3.5.6.2) which is suitable 
for a cross-correlation analysis. 

5.2.1 Scanning the sound wave signals 

Although the illumination and exposure signals must be carefully synchronised 
with each other, the freedom remains to perform the recording at an arbitrary 
delay time after the start of the sound wave signal is sent to the loudspeaker. This 
implies that, by progressively increasing this delay, successive PIV measurements 
can be made to 'scan' the sound wave signal through time and reconstruct it later 
by data analysis. This procedure may be applied at any position along the tube 
to study the spatial variation of the sound wave signal. 

Let t s be the scanning duration of the sound wave under investigation and 
n s be the number of points within the scanning duration that velocity measure- 
ments are performed. Assuming measurements are taken at equal intervals 6t s , 



known as the scanning step, as illustrated in Figure (5.5!, the scanning duration is 
given by 

t s = n s 6t s . (5.2.1) 

The velocimetry time t v , which is the time separation between the start of the 
scanning and the mid-point of the delay between two illumination pulses, may 
be determined as 

t v = u6t s , ^ u < u s , (5-2-. 2) 

where n is an integer. By progressively increasing the quantised time t v , while 
performing a PIV measurement at each time, the velocities associated with the 
sound wave may be obtained at the scanning frequency. 



4 Transistor-Transistor Logic (TTL) signals are commonly used for triggering. The signal is 
in the form of a square-function and in the high mode its value is between +2 V and +5 V. 
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Figure $.$: The sound wave scanning process. 



In order to ensure that all the sound wave information in the scanning dura- 
tion is measured, according to the sampling theorem, the scanning rate u s /t s = 
1/6t s must be at least equal to the Nyquist rate of the sound wave within the 
scanning duration. In other words, if f max is the highest frequency of the scanned 
sound wave, the lower-limit of the scanning rate is 

min(l/&T s ) = 2f max . 

It may be shown that the upper-limit of 1/6t s is determined by the spec- 
ifications of the hardware used for the measurement. Let T; be the minimum 
duration of the illumination pulses which crucially depends on the light source 
illumination power, the recording device sensitivity and the DAQ precision. The 
minimum displacement which can be measured by two identical illumination 
pulses of duration x„ assuming that the displacement is measured between the 
mid-points of the pulses, is within T;. This implies that the input sound wave 
may be measured with the minimum precision of Ti or 

max(l/6T s ) = 1/t ; . 

For the experimental set-up used in this study, the upper-limit of 1/6t s was ob- 
tained to be ~ 10 kHz. 

When carrying out acoustic PIV measurements, the scanning frequency should 
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be set to the maximum value obtainable for the system. Although setting the 
scanning frequency higher than the Nyquist rate does not provide any more in- 
formation about the reconstructed signal, which is assumed to be the same as 
the input signals, by doing so any unexpected behaviour of the output can be 
sampled. 

5.3 Software: Output Data Analysis 



A diagram describing the data analysis procedure is given in Figure ( 5.6 1. The im- 
age evaluation process begins by reading a pair of raw PIV images downloaded 
from camera via the frame-grabber. The images are then pre-processed using 
Gaussian smoothing which is effective in removing background noise. The ker- 
nel size of the smoothing filter was set to 3 x 3. 

Three different algorithms were tried for the main data processing. The FFT- 



based traditional cross-correlation, as described in § 3.5 , was implemented as the 



standard approach. Introduced by Westerweel et al. (1997), discrete window 



offset was used in a second attempt to improve the accuracy in displacement 
detection in the correlation plane. The method is based on spatial offsetting 
between the interrogation windows corresponding to the image pair. The offset- 
ting value is an integer-pixel vector which is chosen to be roughly equal to the 
expected displacement. Hence, the method works best with some pre-existing 
knowledge of the fluid motion, however, its 'built-in' noise removal is beneficial 
in reducing detection error. 

Both of the previous processing methods generate unavoidable erroneous ve- 
locity vectors caused by non-ideal characteristics of the input PIV images. Con- 
sequently, in order to obtain reliable results, applying data validation is neces- 
sary when these methods are in use. The PIV image evaluation algorithm of 
correlation-based correction 'automatically' filters spurious velocity vectors while 
calculating displacements in the interrogation area of interest. The method, 



founded by |Hart (2000), involves using two (or more) interrogation areas for 
detecting each vector. The windows are set apart by a fraction of the interroga- 
tion size. The corresponding correlation plane is yielded by element-by-element 
multiplication of the interrogation areas. 

The post-processing filter was implemented to be used in conjugation with 
data processing methods with inaccurate output. A two-dimensional spatial 
adaptive local median is a common post-processing filter in which each vector 
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is validated by the median of its neighbours. Typical acoustic PIV measurements 
include more than two dimensional data, e.g. two dimensions in space and one 
more dimension in time or sound intensity. Based on this fact, the adaptive local 
median filter was extended to three dimensions using a 3 x 3 x 3 kernel. 
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Figure $.6: Output data analysis diagram. 



Chapter 6 



Measurements 



Because things are 
the way they are, 
things will not stay 
the way they are. 

Bertolt Brecht 



6.1 Sound Pulses 

As a preliminary acoustic PIV measurement, the propagation of a sound pulse in 
tube A was studied. The structure of this experiment is similar to methods for 
finding the impulse response of a loudspeaker with the difference that here the 
propagation of the velocity wave, instead of the pressure wave, is measured. 

6.1. i Generation and Propagation of Sound Pulses 

An ideal pulse may be represented by a 6-function. There are two main factors 
which cause the pulse created in practice to deviate from the ideal pulse: the 
pulse generation device and the medium in which the pulse propagates. 

A wide range of approaches for generating an ideal sound pulse has been re- 
ported in the literature including the use of chemical reactions, spark discharges, 
electro-acoustic transducers and piezoelectric devices which can produce pulses 



of the duration of 0.07 to 2 ms. A review of these means may be found in Mar 



shall 



(1990 1 where producing sound pulses of the duration of 0. 1 ms by a digital 



signal generator, coupled to different sound generation devices, is presented. 
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While the pulse generation device diverges the shape of the created pulse 
temporally, the spatial deviation from the ideal shape is due to the attenuation 
and dispersion properties of the propagation medium. The acoustic harmonic 



motions are known to be decayed spatially in an exponential way (Morse and 



1NGARDII968, §6.4). In contrast, Learned (1979 1 showed that the spatial decay of 



sound pulses is non-exponential and is in the form of power-law ( Kalinichenko 

ET AL.||2QOl[ §5.3). 



6.1.2 Measurement 

As a basic estimation of a 6-function, a square-function of duration of 0.1 ms was 
chosen as the input signal for driving the loudspeaker. Two TTL illumination 
pulses of the duration of 0. 1 ms separated by a delay identical to the illumination 



durations were used to record the PIV images. Figure ( 6.1 ) illustrates the design 
of the input signals. 
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Figure 6.1: Input signals for the measurement of a sound pulse. 



The measurement was executed at a position 18 cm away from the loud- 
speaker and for a scanning duration of 1 ms. At the first attempt, a scanning 
resolution of 0.1 ms was applied to produce 10 PIV images. This makes the 
temporal overlapping between two consecutive scans to be 50%. Spatially, the 
images almost covered the width of the tube. The images were then analysed 
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using correction-based correlation with interrogation areas of size of 256 x 256 



pixels. The results are shown in Figure ( 6.2 ) 



Figure ( 6.3 ) presents the measured velocity pulse obtained by spatial averaging 
of the acquired vector-maps across the tube, assuming that the change of the 
velocity vectors across the tube is negligible. While the reconstructed velocity 
wave is not of a desirable resolution, it confirms that the input square function 
is measured as a decayed bipolar pulse. The bipolarity of the result reflects the 
response of the loudspeaker to the input square function. 

In order to have a better picture of the pulse, as well as testing the system at 
a high scanning rate, a second measurement was carried out with an improved 
temporal resolution of 0.01 ms, keeping the same scanning duration of 10 ms. 
Hence, the temporal overlapping of two consecutive scans in this measurement 
is 99%. This produced 100 PIV images which were analysed with the same 



parameters as the previous measurement. The result is given in Figure (6.4 1. 

It is worthy of notice that the extremely high temporal overlapping of the sec- 
ond measurement not only did not generate over-sampled repeated information, 
but dramatically improved the reconstructed velocity pulse. Using this data, an 
excellent fit in the form of a non-symmetrical bipolar Gaussian was obtained 



describing the impulse response of the loudspeaker (Figure 6.4). Figure (6. 5 1 
compares the Fourier transforms of the measured velocity pulses with different 
scanning resolutions. 



6.2 Linear Harmonic Oscillation 

The acoustic velocimetry of standing sound waves inside the tube was carried 
out with the tube driven at its first natural frequency. The tube was coupled by a 
loudspeaker at one end and rigidly terminated at the other end (§ J5.i.i| ). Using an 
oscilloscope, the first two resonance frequencies were determined to be around 
280 Hz and 570 Hz corresponding to periods of about 3.6 ms and 1.8 ms. 



Figure (6.6 1 shows the input signals used for this measurement. A sinusoidal 
sound wave, with a period matching the resonance frequency of interest, was sent 
to the loudspeaker while the illumination signal scanned the sound wave over a 
full acoustic cycle using a temporal scanning resolution of 0.1 ms. This scanning 
resolution yields a sampling frequency of 10 kHz which is around 18 and 9 times 
the Nyquist frequency of the two first modes, respectively. 

89 



Measurements 



t v = 0.1 ms 



t v = 0.2 ms 
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t v = 0.4 i 
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t„ = 0.6 ms 
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10 15 20 25 30 
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t v = 0.8 ms 
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Figure 6.x: Velocity vector-maps illustrating the temporal scanning of a sound pulse 
propagating along the tube. The measurements were carried out 18 cm away from the 
loudspeaker. The correlation-based correction algorithm was used with interrogation 
areas of size of 256 x 256 with 93.75% overlapping. The horizontal axes show the 
width of the tube in mm. 
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Figure 6.3: Transversal mean acoustic velocity variation of a sound pulse propagating 
along the tube with the scanning resolution ofO. 1 ms. Each point corresponds to the 
mean value of the velocity vectors in the vector-maps of Figure (6.2). 



3 

•;■: 




Figure 6.4: Transversal mean acoustic velocity variation of a sound pulse propa- 
gating along the tube with the scanning resolution of 0.01 ms. The solid line rep- 
resents a bipolar Gaussian fit as a function of time given by i.29e~ (0 - 87t ~ 2 - 86) — 
0.34 e -(°- 35t - 2 - 56 ) 2 . 
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Figure 6.$: Comparison of the Fourier transforms of the velocity pulse with scanning 
resolutions of 0.1 ms (10 samplings, dashed line) and 0.01 ms (100 samplings, solid 
line). 



The sound pressure level at the rigid termination was measured by a probe 
microphone (Bruel & Kjasr, 4170) connected to a measuring amplifier (Bruel & 
Kjasr, 2607). The needle of the probe microphone was entered into the tube 
through a small hole which was made at the end of the tube. Since the input 
sound intensity was controlled by driving the loudspeaker using a Hi-Fi amplifier, 
the probe microphone was applied to 'calibrate' the Hi-Fi amplifier. Results are 
given in Table ( |6.i| . 

6.2.1 Variation with Time 

Driving the loudspeaker at the two first natural frequencies of the tube, the sound 
pressure level at the rigid end of the tube was read from the measuring amplifier 
to be 124.0 dB for the first mode and 123.4 dB for the second mode. The 



SPL values are corresponded to the Hi-Fi amplifier being set to 2x (Table 6.1 1. 
Monitoring the acoustic pressure on an oscilloscope using the probe microphone 
confirmed sinusoidal variation of acoustic pressure. 

Based on the described input signals, the PIV measurements were performed 
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Sound Wave 
Illumination:Begin 

Illumination:End 
-2 ( 




Time (ms) 



Figure 6.6: Input signals for the measurement of standing waves. 



Hi-Fi amplifier dial SPL of the first mode (dB) SPL of the first mode (dB) 



2x 142.0 
4x 151.2 
6x 155.0 



143.4 
153.0 
156.9 



Table 6.i: Calibration of the Hi-Fi amplifier using the probe microphone. The SPL 
values are measured at the rigid end of the tube. 



at a distance 1 5 cm away from the loudspeaker. The PIV images were analysed 
by interrogation areas of size of 256 x 256 using the correlation-based correc- 



tion algorithm. This algorithm generated vector-maps, shown in Figure (6.7 1 



which did not need any post-processing including data-validation for excluding 
the spurious vectors. 

Averaging the resulting acoustic velocity vectors across the width of the tube, 
the temporal variation of the standing wave was obtained as shown in Figures 



(6.8a) and (6.8b) for the first and the second modes, correspondingly. Both 
variations are almost in a sinusoidal form with an acoustic velocity amplitude of 
about 0.66 m/s (first mode) and 0.48 m/s (second mode). As the experiment 
was not designed to extract any temporal phase information, throughout, in 
presenting the results the temporal phase is manually set to zero. 



Figure (6.8c I shows the Fourier components of the measured acoustic veloc- 



ity for the first two modes. The two highest peaks in the frequency domain are 
around the two natural frequencies of the tube. The contribution of the rest of 
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t T = 1.6 i 



t T = 1.7 ms 



t v = 1.8 ms 
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Figure 6.7: Velocity vector-maps of the harmonic oscillation at the first resonance 
of the tube with a SPL of 126.2 dB and at a position 15 cm away from the 
loudspeaker.^ 
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t„ = 1.9 ms 



t v = 2. ms 



1/niNiNiNiniiii 



t v = 2.1 ms 



immiiiiini 



2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 
t v = 2.2 ms t v = 2.3 ms t v = 2.4 ms 



, 



2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 
t v = 2.5 ms t v = 2.6 ms t v = 2.7 ms 







2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 
t v = 2.8 ms t v = 2.9 ms t v = 3. ms 
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t v = 3.1 ms t v = 3.2 ms t v = 3.3 ms 

1 1 I I I III III 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11 1 
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t v = 3.4 ms t v = 3.5 ms t v = ms 




2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 2.5 5 7.5 10 12.5 15 17.5 20 



*" The correlation-based correction algorithm was used with interrogation areas of 
size of 256 x 256 with 93.75 % overlapping. The horizontal axes show the width of 
the tube in mm. 



95 



Measurements 



the peaks, especially for the first mode, is quite low. 

6.2.2 Variation with Sound Intensity 

In the previous measurement, the standing wave of acoustic velocity was mea- 
sured for constant values of sound intensity and axial distance. The measure- 
ments were extended by changing the sound intensity while keeping the rest of 
the parameters the same. 

At 23 cm away from the loudspeaker (0.19 Ai, 0.37 A2), the acoustic velocity 
inside the tube was temporally scanned for the first two modes when the Hi-Fi 
amplifier was set to 2 x, 4 x and 6x (Table pO). Analysing the PIV data, based 
on the same parameters used previously, the transversal mean acoustic velocity 



was obtained as presented in Figures (6.9a) and (6.9b! for the two first modes 



respectively. It is noticeable from the results that as the sound intensity increases 
inside the tube, the measurement system fails to gain expected values. Later, in 
this work, the sources of this anomaly will be investigated, however-for the time 
being-the results of this investigation is used to analyse the measurements. 

The algorithm employed so far in the processing of the PIV images is founded 
upon Correlation-Based Correction (CBC). It will be shown that using correlation- 
based correction together with the Symmetric Phase-Only Filter (CBC-SPOF) 
dramatically improves the anomalous results towards their expected values. 



Figures (6.9c) and (6.9dl demonstrate the variation of acoustic velocities by 
processing the same PIV images using CBC-SPOF, which imply a significant 
change when compared to the results obtained by CBC processing. 

Calculating their RMS values, the variation of the measured acoustic veloci- 
ties at the given axial position (0.19Ai, O.37A2) versus the measured SPL at the 



rigid end is shown in Figure (6.10 1. As illustrated, the CBC-SPOF processed 



results well-agree with the theoretical prediction. 

6.2.3 Variation with Space 

Extending the measurements in one more dimension, the former experiment was 
carried out along the axis of the tube. In detail, the standing acoustic velocity 
wave was temporally scanned for different frequencies, sound intensities and po- 
sitions forming the standing wave patterns of the first two modes. 



Figure (6. 11 1 shows the results analysed by CBC-SPOF processing. A com- 



parison of the experimental values with the theoretical standing wave pattern for 
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a) Acoustic Velocity of the first mode 
with the SPL of 124.0 dB. Each 
point corresponds to the mean value 
of the velocity vectors in the vector- 
maps of Figure \6.j). 



b) Acoustic Velocity of the second 
mode with the SPL of 123 .4 dB. 










c) The Fourier transforms of the 
acoustic velocities of the first mode 
(solid line) and the second mode 
(dashed line) of the tube. 



Trequency (kHz) 

Figure 6.8: Transversal mean acoustic velocity variation of the first two modes of the 
tube and their Fourier components. Measured at 15 cm away from the loudspeaker 
with the scanning resolution ofO. 1 ms. 
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6.2. Linear Harmonic Oscillation 
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Figure 6.10: Variation of RMS acoustic velocity at 25 cm away from the 
loudspeaker with SPL measured at the rigid end. A comparison between 
the data analysed by CBC (f) and CBC-SPOF (k) is given together with 
the theoretical values ( — ■). The RMS acoustic velocities are calculated by 
the values presented in Figure (6.9). 
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Figure 6.11: Variation of RMS acoustic velocity as a function of the axial 
position. All values are obtained by CBC-SPOF processing. The solid 
lines show the corresponding theoretical standing wave pattern of a tube 
with perfectly rigid terminations. The loudspeaker is located atx = cm 
and the rigid end is atx = 62 cm. 
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6.3- Streaming 



a tube with perfectly rigid terminations is also given. The results around the 
anti-node are in a good agreement with the predicted values. However, the mea- 
surements around the nodes are biased towards lower velocities. It seems that 
this bias is independent of the absolute value of the acoustic velocity. 

Another type of deviation from the theoretical values is evident at positions 
close to the loudspeaker (x = cm) where there is an acoustic velocity drop. This 
may be explained by pointing out the fact that as the loudspeaker termination is 
not absolutely rigid, there is no acoustic velocity node around this position. This 
implies that the standing wave pattern around the loudspeaker has a distorted 
form. 



6.3 Streaming 

6.3.1 Steady Streaming Visualisation 

In order to investigate the presence of streaming, a measurement was carried out 
to visualise the tracer particles on a time-averaged basis. 

Unlike the other measurements in this work, this experiment is a single- 
frame multiple-exposure measurement. An exposure signal of the length of a 
few acoustic cycles is sent to the camera (where each acoustic cycle is determined 
by the frequency of the sinusoidal sound wave signal). During the exposure sig- 
nal, illumination pulses at the frequency of the sound wave signal are sent to the 
laser. For a better resolution, a lens of higher magnification was used (Lens B in 



Table (5.1!). The result is a super-imposed collection of particle-images. 

In order to verify the reliability of the measurement, this test was carried out 
on an object connected to an electrical shaker vibrating at a frequency of 2 kHz. 



Figure (6.12a) shows the result for the first 20 cycles of the vibration of the object 
when it starts to oscillate. It may be seen that after around 10 cycles the object 
reaches equilibrium where its image intensity is higher. 

A similar procedure was applied to the tracer particles in the first resonance 
of the tube around the anomalous location along the tube. In order to ensure 
that the sound field was in its equilibrium state, the measurement was taken after 
100 acoustic cycles. To avoid any streaks, each illumination signal was synchro- 
nised to be sent around the phase position corresponding to zero acoustic veloc- 
ities. Figure (6.i2b| shows the super-positioning of the particle-images for 10 



acoustic cycles. The expected super-position pattern would be a set of regularly- 
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6.4- In-Plane Noise 



spaced particle-images corresponding to the number of recorded acoustic cycles. 
However, unexpectedly, this number for each individual particle-image does not 



exceed more than 2. The cause of this result will be investigated later in §6.5.2.2 



6.4 In-Plane Noise 

An important outcome of the previously presented acoustic PIV measurements is 
the deviation of the experimental acoustic velocities from their expected values. 
The deviation is observed to be a function of space, time, and sound intensity but 
not the predicted or measured value of the acoustic velocity. The latter implies 
that the deviation cannot be immediately linked to the upper-limit restriction 
of the acoustic PIV system. In order to explain the anomalous behaviour of the 
measurement system, two main questions are investigated. First, an attempt is 
made to address the source of the deviation. Secondly, based on the cause of the 
anomaly, means to minimise the deviation are demonstrated or suggested. 

6.4.1 Probing Correlation Plane 

Since the information associated with the motion of the fluid under a PIV mea- 
surement is stored in the correlation plane, it is sensible to start tracing the source 
of the deviation by examining the correlation plane. 

Figure ( 6.13) ) certifies the correlation planes obtained by some acoustic PIV 



measurements presented previously, namely, a sound pulse and a standing wave 
with different values of sound intensities. One common feature shared between 
the demonstrated correlation planes is the broad peak along the direction of 
illumination. It will be shown that the presence of this wide peak is due to 
the narrow height of the light sheet. 



The effective height of the light sheet, measured by eye as given in §5.1.3 



was around 1 cm. Figure (6.14)) shows the illumination profile, visualised by the 



seeding particles, as recorded by the CCD camera. The illumination duration 
of this recording was 0. 1 ms which is the same duration used in the presented 
measurements. The image indicates that the effective height of the recorded 
pulsed illumination profile is of the order of 1 mm which is about 10% of the 
one of continuous illumination profile as measured by eye. One reason behind 
this is the lack of sensitivity of the CCD chip to low-power light. As a trade- 
off, this deficiency may be compensated by increasing the illumination power. 
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a) Correlation plane of a sound 



pulse (Figure 6.4), 



b) Correlation plane of the ex- 
tremum of a standing wave (Figure 
'^ with a SPL of 142.0 dB. 




c) Correlation plane of the ex 
tremum of a standing wave (Figur 
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d) Correlation plane of the ex 
tremum of a standing wave (Figur 
Tyc} with a SPL of 1 55. .0 dB. 
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Figure 6.13: Comparison between the experimental correlation planes 
corresponding to acoustic motions. All correlation planes are obtained by 
CBC processing. The brighter regions stand for higher values. 



6.4- In-Plane Noise 




Tube Length (mm) 

Figure 6.14: The illumination profile as recorded by the CCD camera with an illu- 
mination duration ofiO. 1 ms. 



However, the laser used in this work had an output power of 250 mW which is 
classed as low to moderate power. Typically, lasers applied in PIV measurements 
are pulsed lasers with an output energy of 50 to 120 mj per pulses of duration 
of 5 ns which is associated with a power of about 10 MW. Accordingly, while it 
is not possible to directly compare the output power of a continuous laser to a 
pulsed one, it is reasonable to assume that the laser used in this study is of many 
orders less powerful than a recommended laser for PIV measurements. However, 
the fact that the developed set-up, based on the low-power laser, was able to 
measure velocities up to a few meters per second is significant. 

The narrow height of the recorded light sheet, relative to the size of the 
interrogation areas, causes a strong illumination gradient across the interroga- 
tion areas. Consequently, this induces a background correlation peak H^ c { s ) 



(see §3.5.6.2) to appear as broad peak. The background correlation peak has a 
Gaussian form with a width of l/8th of the interrogation area size. 



6.4.2 Filtering Correlation Noise 

The discussed illumination and recording limitations do not prevent the data 
acquisition system from losing the displacement signal. In fact, the displace- 
ment peak signal is buried under the background correlation noise and may be 
extracted by an appropriate filter. 

Removing the noise of correlation peaks is neither a new idea nor restricted to 
PIV data processing. Optimising the correlation signal-to-noise ratio was stud- 
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ied in the 1970s where it was applied to regain radar signals from noisy data 
( Carter|i976 ; Knapp and Carter|i976| ). The approach was first applied to PIV 
correlation planes by Wernet and Edwards (1988) based on optical correlators. 
More recently, an extension is implemented by Wernet ( 2005I where the author 
demonstrates the effectiveness of phase-only filters (POF), in particular Symmet- 
ric phase-only filter (SPOF), on PIV image evaluation. The basis of phase-only 
correlation is matching by the phase of the Fourier transformed data. This is in 
contrast to the traditional method which tries to match based upon the ampli- 
tude of the frequency domain. The significant advantage of phase matching over 
amplitude matching is that the phase component stores information about the 
position of the input objects, while the amplitude component saves information 
about their size. Correlation obtained by position-only matching not only yields 
a stronger displacement peak, but also fades the rest of the peaks remarkably, 
which may be counted as noise. Since the signal of interest in a PIV analysis is 
the displacement peak, SPOF effectively eliminates the unwanted peaks includ- 
ing background correlation. 

The application of SPOF in recovering acoustic displacements measured by 



low-light PIV is manifested in this work (§6.2 1. Figure (6.15a! illustrates a cor 



relation plane obtained by CBC which uses Fourier amplitude matching. The 
correlation plane, corresponding to a PIV image recorded by the previously men- 
tioned narrow illumination, clearly contains a wide background correlation peak 
which conceals the displacement peak. The displacement peak, as shown in Fig- 
ure ( 6.15b)), is retrieved by using phase-only correlation. 



6.5 Out-of- Plane Noise 



6.5.1 Acoustic Motion Visualisation 

A sound field visualisation experiment with high-magnification was set up to 
record the microscopic motion of the particle-sources. The aim of the test was 
to provide a direct observation of the motion of the particle-sources under the 
anomalous condition using the existing PIV set-up. The specifications of the 



measurements are as described in §6.2.1 the only difference being the use of a 



camera lens of greater magnification (Lens B in Table (5.1!) 



Figures (6.16! and (6.17) show visualisations of the sound field generated 
when the tube was driven at its first natural frequency. The visualisations were 
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Sx 



a) CBC 




Sx 



b) CBC-SPOF 



Figure 6.15: Comparison between correlation planes obtained by Fourier amplitude 
matching (CBC) and phase-only matching (CBC-SPOF). The three-dimensional 
plots correspond to the indicated square regions about the displacement peaks. The 



correlation plane is the same as in Figure (6.13 b) 
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Figure 6.16: Visualisation of sound field in the non-anomalous location 
108 °f the tube. 
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Figure 6.17: Visualisation of sound field in the anomalous location of the 
tube. 
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carried out at, respectively, the non-anomalous and anomalous positions inside 
the tube and at selected times within an acoustic cycle. For ease of observation, 
negatives of the images are presented. The visualised sound field is comparable 
to that presented in the works of Brandt and Hiedemann (1936)) and Hann 



(199 5 1 (see Figures ( 4.3a! and (4.3b) on page 61); however, the temporal resolu- 
tion achieved in this study is of the order of 1 00 times greater than that of the 
mentioned works. Consequently, instead of having a 'dumbbell-shape' recorded 
over a few acoustic cycles, each particle-image is expected to appear as a streak 
at the extremum points of the acoustic velocity and as a 'half-dumbbell-shape' as 
they get closer to smaller velocities. 

A qualitative comparison between the two sets of visualisations shows that 
the streaks recorded under the anomalous condition are relatively shorter and of 
stronger intensities than the non-anomalous streaks corresponding to the same 
phase. This is an unexpected result as the anomalous streaks, whether associated 
with higher sound intensities or regions around the velocity node, should be 
generally longer than the non-anomalous streaks. A more surprising outcome is 
observed about the temporal extremum points, t = 0.3/f and t = 0.8/f . At these 
points the non-anomalous streaks have a predictably uniform intensity profile. 
However, the anomalous streaks of the same phases generated an intensity profile 
which looks like full-GAussiAN or partial-GAussiAN. 

6.5.2 Discussions 

The unexpected behaviour of the visualised acoustic motion may be interpreted 
once assuming that, during the illumination time, the particle-sources were laid 
on a plane parallel to and within the light sheet volume; and once assuming that 
they have a considerable displacement normal to the light sheet while they were 
recorded. 

6.5.2.1 In-Plane Motion 

The first of the above assumptions implies that the recorded streaks reflect the 
true two-dimensional motion of the particle-sources. Accordingly, it may be 
argued that particle-sources actually do not perform the expected harmonic os- 
cillations. One may go further and expand this possibility to whether the source 
of the unpredicted motion is from the fluid flow itself or due to the fact that the 
particle-sources do not follow the flow faithfully. 
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The former hypothesis was rejected by measuring the acoustic pressure using 
a probe microphone. The measurement, which was carried out similarly to the 



description in § 6.x confirmed the sinusoidal form of acoustic pressure variation 
implying the same result for acoustic velocity. 

Under the measurement conditions, given a frequency of the order of 10 2 Hz 
and SPL of less than about 160 dB, the infidelity of the particle-sources to the 
fluid flow in the sound field does not agree with the classical theory presented 



in §4.5.2.2 Even by neglecting the measurement conditions, the classical theory 
simply predicts the infidelity in the sense of attenuation in the amplitude of 
oscillation. This makes tracer particles infidelity to be a less-likely source of 
anomaly. 

6.5.2.2 Out-of-Plane Motion 

The idea of the motion of the particle-sources during the illumination time may 
be supported by a number of reasons. First, the stronger intensity of the anom- 
alous particle-images may imply their motion normal to the recording plane, 
therefore, creating shorter streaks. Secondly, the unexpected Gaussian pro- 
files of the anomalous streaks, at the temporal extremum points, can prove that 
the particle-sources cross the light sheet volume. Consequently, the Gaussian 
profile of the light sheet thickness is 'mapped' to the expected uniform profile 
of the streaks making them appear as Gaussian streaks. Hence, according to 
whether the particle-sources cross the light sheet volume completely or partially, 
they generate streaks with full-GAussiAN or partial-GAussiAN profiles, respec- 
tively. Finally, considering the out-of-plane motion of source particles explains 
the streaming visualisation pattern discussed formerly (Figure 6.i2a| by the fact 



that the time-averaged motion normal to the light sheet cuts the continuation of 
the regularly-spaced particle-images. 

One reason behind the out-of-plane motion could be due to the angle be- 
tween the light sheet normal and the tube axis. In the PIV set-up used for the 
measurement of acoustic oscillatory motions, it is assumed that the normal vec- 
tor of the light sheet is perpendicular to the axis of the tube. This prevents all the 
particle-sources initially located within the illumination volume from exiting the 
volume while avoiding new particle-sources simultaneously entering. However, 
in practice this angle may deviate from 90°. 



Let the deviation angle be denoted by d. From Equation (3.5.25), for a 
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Figure 6.18: Variation of out-of-plane loss-of-correlation with the deviation of the 
angle between the Gaussian light sheet normal and the tube axis for different values 
of thickness of light sheet. The acoustic motion corresponds to an extremum point in 
space-time with v = 2 m/s and%i = 0.2 msfor the first mode of the tube. 
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Figure 6.19: Variation of out-of plane loss-of-correlation with time atd = 5 °, vo 
2 m/s andxi = 0.2 msfor the first mode of the tube. 
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Figure 6.20: Variation of out-of-plane loss-of-correlation with SPL at the velocity 
node of the first mode of the tube with d = 5 °andii = 0.2 ms. 
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Figure 6.21: Variation of out-of plane loss-of-correlation with axial position at a 
temporal extremum with d = 5 °and%i = 0.2 ms. 
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Gaussian intensity profile of the light sheet, the out-of-plane loss-of-correlation 
F may be given as 

V 71 a / o f x + ( v ° sin(k n x) sin(cu n t))Ti sin(-&) N " 



F (Az L ) = ^Az L expl-8^^^ v - '^ n " l ^j j , (6.5.1) 

where v is the magnitude of the acoustic velocity, T ; is the duration of the il- 
lumination and Azl is the thickness of the light sheet. Since the determining 
experimental parameter in this equation is the light sheet thickness, it is worthy 
to analyse the behaviour of F Q as a function of Azl. 



Figures (6.18 1 to (6.21) examine the variations of F (Azl) as a function of 
deviation angle, time, SPL, and axial position, respectively. In order to apply 
these results to the previous measurements, Fo is calculated for the first mode 
of the tube under investigation with an illumination duration of 0.2 ms. The 
graphs show a significant drop of valid detection of the correlation peak for an 
effective light sheet thickness of about 0. 1 mm and less. Accordingly, the valid 
detection probability decreases at some points in space-time-intensity, namely 
around the extremum points of the phase and the velocity node of the tube. 
This is in general agreement with the previous unexpected observations. 

Based on the calculated out-of-plane loss-of-correlation, and assuming that 
the in-plane loss is negligible, the anomalous measured velocity may be described 
as 

v* (x,t,Az L ) = F o (Az L )v sin(k n x) sin(u> n t). (6.5.2) 



Figure (6.22) compares this theory to the experimental data presented in §6.2.1 
The fitted curve is associated with a Gaussian light sheet with a thickness of 
about 0.05 mm and deviation angle of 9°. Since the height of the light sheet, 
as recorded by the camera, was observed to be around 10% of its actual value 



(§6.4.1), it is reasonable to assume that the same ratio is applicable to the thick- 
ness of the light sheet. 

It is possible to decrease or eliminate the out-of-plane noise using a more 
powerful light source. In this case, the main benefit of more light energy is 
to produce a thicker light sheet and not obtain higher quality particle-images. 
Another option, which is harder to achieve, is to carefully align the laser sheet 
relative to the axis of the tube. Moreover, the use of a camera with a higher 
sensitivity allows shorter illumination pulses to be achieved without a significant 
loss of the effective light sheet thickness. 
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Figure 6.22: Comparison between CBC-processed experimental acoustic velocities 
(•), the theoretical fit considering the out-of-plane effect ( — ) and the expected values 
( — ) for the first mode of the tube with SPL = 151.2 dB and T; = 0.2 ms. The 
fitted curve corresponds to a Gaussian light sheet with a thickness of 0.05 mm and a 
deviation angle of 9 °. 
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7.1 Summary and Conclusions 

7.1. 1 Introducing acoustic PIV 

While a lot of research has been carried out on the acoustic velocimetry of low- 
speed flows, there is a perceptible lack of literature in measuring acoustic ve- 



locities particularly by non-intrusive whole-field means (§4.3!. This study is an 
attempt to fill this absence by introducing the technique of acoustic PIV. 



The verification of the principles of PIV (§3.2.3) in the presence of a sound 



field has been carried out by examining various interactions between the mea- 



surement components and the sound field (§4.4). Among these interactions, the 
effect of sound fields on the tracer particles was found to be negligible for low to 
moderate acoustic frequencies and sound intensities. 



7.1.2 Scanning Sound Waves 

The main idea behind the acoustic PIV experimental set-up is to be able to mea- 
sure the velocity of the sound wave of interest at a high-enough temporal rate. 
An absolute scanning rate of up to 5 kHz has been achieved in this study, which 
is suitable for the precise measure of the time variation of sound waves within an 



acoustic cycle (§6.2.11. 

The sound wave scanning process has been successfully applied to both sound 



pulses (§6.1 1 and sinusoidal waves (§6.2). Variation of acoustic velocities with 
sound intensity and the standing wave pattern of the first two modes of a simple 
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tube has been studied. In addition, more improved visualisations of acoustic and 
hydrodynamic motions were presented. 

7.1.3 Low-Light High-Speed PIV 

The laser used in this study was a CW solid-state laser externally triggered to pro- 
duce quasi-pulses. The power of this laser was of many orders lower than typical 
pulsed lasers used in PIV measurements. Considering the fact that measuring 
acoustic velocities with a standard PIV system is near to the limitation border of 
the technique, it is a challenging task to perform acoustic PIV with low-power 
light sources. 

Two main problems were encountered using the low-light system. The first 
one was due to the narrow height of the Gaussian light sheet used in this study 
which caused a significantly broad background peak in correlation planes. This 
noise, referred to as the in-plane noise, was successfully eliminated using a phase- 



only PIV image processing method (§ 6.4). The second problem was the presence 



of out-of-plane noise which was caused by the narrow thickness of the light sheet. 
It was demonstrated that this noise led to a excessive loss of valid detection at 



specific space-time-intensity points (§6.51 



7.2 Future Research 

7.2.1 Instantaneous Measurement of Impedance 

One immediate application of acoustic PIV is the calculation of acoustic im- 
pedance. While it is possible to measure the acoustic pressure independently 
using conventional methods, it is desirable and more convenient to calculate the 



pressure from the velocity vector-map (Imaichi and Ohmi 1983). Combining 



this method with acoustic PIV has the advantage of measuring variation of im- 
pedance simultaneously and non-intrusively Estimating the acoustic pressure via 
the acoustic velocity is based on the numerical calculation of the spatial gradient 
of the velocity, implying that the velocity vector-map must be whole-field with 
valid representations of the sound field velocity values in each interrogation area. 
Evidently, this requires the tracer particles to be in the absolute fidelity phase. 

Measurement of the impedance in this way has potential applications in noise 
reduction or the study of musical instruments. One attractive benefit of employ- 
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ing acoustic PIV in these applications is that the technique is independent of the 
shape of the flow passage, enabling acoustic measurements to be carried out in 
complex geometries. 

7.2.2 Aerosol Measurements 

Although this work is based on using aerosol particles to measure sound fields, 
acoustic PIV may be employed to use sound fields in the measurement of aerosols. 
Velocity vector-maps corresponding to sound fields store more information than 
just the acoustical properties. Indeed, the velocity distribution extracted from an 
acoustic PIV measurement contains data about the size distribution of the tracer 



particles (§4.5.3.21. In addition, it is possible to study the motion and internal 



interactions of aerosol particles in sound fields using acoustic PIV. 



119 



Ill 



Appendices 



121 



Appendix A 



Time Series 



Broadly speaking, physical phenomena are processes varying in time, e.g. 
'wave-like' structures. These processes are known as time series. In the 
following appendix, those elements of time series theory which are re- 
ferred to in this study are presented in a compact form. The concepts are mainly 



taken from Priestly (1981) 



A.i Terminology 

Definition A. i. Any functioning which can yield one of many likelihoods is 
called an experiment <£. 

Definition A.ii. All possible outcomes of an experiment (S are referred as the 
sample space O of that experiment. Each single possibility in a sample space is 
named an element CI. 

Definition A.iii. Any sub-set of a sample space O corresponding to an experi- 
ment <£ is defined as an event £ of that experiment. 

Definition A. iv. Assume that there is a mapping from every event £, belonging 
to a sample space D., to a real number p (£). The number p (£) is the probability 
of the event £ and fulfills the following conditions 

0^p(£Kl (Ai.ia) 



p[Cl) = l; (Ai.ib) 
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Definition A. v. For every element Cl, belonging to a sample space Q., there exists 
a unique real number X[C1) called a random variable. 



A. 2 Univariate Distributions 

Consider, for any real number ?c , the event £(?c ) is given by 

EU)= {CI; (OeQ)A (X(C1) ^ * )}. (A.2.1) 

The univariate distribution function f[?c) of the random variable X(C1) is defined 
as the probability of the event £{?t ), i.e. 

FU)=p[EU)]. (A.2.2a) 

The distribution function may also be expressed explicitly in terms of the corre- 
sponding random variable: 

?U)=p[X^x.], (A.2.2b) 

where p [X ^ %] is the probability for the random variable X. 

A.2.1 Discrete Distributions 

Assume X is a discrete random variable whose possible values are restricted to the 
set {?c j, ?c 2 > K. 3> ■ ■ •}• The discrete probability p t for X has the following properties 

p.=p[X = ? Ci \, (A.2.3a) 

0<p t <l, (A.2. 3 b) 

X?i =1 ' ( A - 2 -3 c ) 
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A. 2. Univariate Distributions 
where i is a natural number. The discrete distribution function for p is given by 

The mean and variance of the discrete distribution function ?{?t) can be calcu- 
lated from 

\l=£_X.iPi> (A.2.5) 



anc 



(T 2 = ^( ?i -U) : 



(A.2.6) 



A. 2. 2 Continuous Distributions 



The random variable X is said to be continuous if we can define a function /(^ 
with the following specifications for all ?c 









(A.2.7a) 



/U)^o, 



(A.2. 7 b) 



+ 00 



fti)dx. =1- 



(A.2.7C) 



/(^) is called the probability density function (PDF of X. Given the PDF of 
a random variable, the corresponding continuous distribution function may be 
calculated by 



7{iL. 



/k)<k. 



(A.2.8) 



The mean and variance of the continuous distribution function !F{?t) is given by 



* + oo 



H 



?cf[?i)d?c : 



(A2.9) 
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and 

a 2 



(x-v) 2 fMd^. (A.2.10) 



A. 2. 3 The Expectation Value 

The statistical properties of a distribution function can be presented in a compact 
form by the expectation operator £{ }. The expectation value of a given function 
g[X) is defined as 

■E{g{X)}=^gU i ) P ., (A.2.11) 



when X is discrete and 



£{g(*)} 



g(a:)/(a:)da:, (A.2.12) 



when X is continuous. 

Using the expectation operator, the definitions of mean and variance can be 
re-written as 

u = -E{X}, (A.2.13) 

and 

a 2 = £{(X-u) 2 }. (A.2.14) 



A.2.4 Special Distributions 
A.2.4.1 Binomial Distribution 

The binomial distribution with parameter n is defined as 



The distribution describes a sample space with only two elements, one of which 
has a probability p of occuring. The discrete probability p i . n represents the 
likelihood of n outcome with probability p . 
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A.3. Multivariate Distributions 



A.2.4.2 Poisson Distribution 



Consider the binomial distribution Equation (A.2.15) when n — > 00 and p — > 
with the condition mp = a, where a is a positive constant. It can be shown that 
p i , with the given assumptions, can be written 

a 1 
?i. a = -e a , i= 1,2,3,..., n. (A.2.16) 



1 



This discrete probability is known as the Poisson distribution with parameter a. 

A.2.4.3 Gaussian Distribution 

Assuming that ^ is a continuous random variable, the associated PDF 

1 U-n) 2 

/k)= ... e ^ 2 , (A2.17) 

Z7l(J z 2 



2 



is said to have a Gaussian distribution function with parameters u ^W a 



A.3 Multivariate Distributions 

The rationale behind univariate distributions can be extended to multivariate 
distributions when describing more than a single random variable. 

A multivariate discrete probability distribution for n discrete values is defined 
as 

PuA 2 ,...,i n =P^ =X. il ,X2 = ?Ci 2 ,...,X n = XLiJ, (A3.1) 

with 

ii = 1,2,... 

12 = 1,2,... 

tn — 1, X, . . , . 

For n continuous variables, the multivariate probability density function is given 
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by 



lim 



flX 1' X 2 -> ■ ■ ■ ' Xn) ~~ 

6*1 ••• ^ n 



(A.3-2) 



Based on Equation (A.2.2D), the definition of the distribution function can 
be extended to its cumulative form: 

?{X\,X2'--->X n ) = P [X\ ^ X\,Xi ^ X 2 '---> x n < * J; (A.3.3) 

where for the discrete and continuous random variable it is given by 

?{Xl,X 2 ,---,X n ) = X Y- ••• X fiLla-ln' ^ A - 34 ^ 



and 

slXiiX2->---->Xi 

respectively. 



— 00 « —00 



f[Xi,X 2 ,---,X n )dXidx 2 ---dx n , 

(A.3.5) 



A.4 Stochastic Processes 

Assume that the outcome of an experiment £ is a quantity X which is a contin- 



uous function of time t. According to definition Equation (A.vl, the parameter 
X[t) is a random variable of the sample space O of <£. Let O have n elements 

{Cl u Cl 2 ,...,Cl n }, 

for each time t the experiment <£ can have a set of 'possible outcomes', called 
records, shown by 

{l(t ) (l 1 ),l(t,Q 2 ),... ) l(t,Q n )}. 

This set of the likely records is named the ensemble. Each single record in the 
ensemble is called to be a realisation and is typically denoted by X{t, D.). 
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A. 4 . Stochastic Processes 

Assume that t varies over all its possible values and for each value a record is 
obtained. In this way, a set of realisations is generated which is referred to as a 
stochastic process 1 . A stochastic process with n random variables is shown by 

{X[t)} = {X{t x ),X{t 2 ),...,X{t n ]\ 

and its corresponding cumulative distribution function is denoted as 

!Fx(t 1 ),X(t 2 ),...,X(t n )l?C 1> ^2>- ■ ■ '^rJ- 

The preceding definitions are also valid for discrete random variables. How- 
ever, in most practical cases, the random variables of a stochastic process is a 
continuous parameter. This will be assumed in the forthcoming sections sec- 
tions. 

It is worth noting that the parameter t in a stochastic process {^"(t)} need not 
necessarily represent time. Depending on the physical implementation of the 
process, it can be any quantity. 

A.4.1 Stationary Univariate Processes 

Let a stochastic process for any ti, t 2 , . . . , t n be written 

{x[u),x[t2),...,x{t n )}, 

and the same process observed with any arbitrary delay or lag t be written 

{JC(ti + x),X[t 2 + t), . . . ,X[t n + t)}. 
The process is said to be stationary if 

9 r x(t 1 ),X(t 2 ),...,X{t n )(?C l»^2>- ■ ■ ' <tn.J = Fx(t l +T),X{t 2 + T),...,X{t n + T)(?C 1>^2> • • • > X. n) ■ 

(A.4.1) 

A stationary stochastic process describes a physical system in which the statistical 
properties do not change over time, i.e. the system is in steady-state or equilib- 
rium. 

For the rest of this study we assume the processes are stationary. 



'The so-called random process. 
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A.4.1.1 Mean and Variance 



Consider the PDF / t (a: ) of the random variable {JK(t)} of a stationary process. 



Using Equations (A.2.81, (A.2.9), (A.2.13) and (A.2.14I, the mean and variance of 



this random variable are given by 

mean{X(t)} = "£{X{t)} 
and 



t/tk)<k = n(t) 



(A.4.2) 



j r + 00 n 

var[JC(t)}=E{(jC(t)-^t)) }= (*-u(t)) / t (* ) d^ = a 2 (t) 

J —OO 



(A.4.3) 



A.4.1.2 Covariance and Correlation Coefficient 

The covariance of two random variables X[t{] and X{t 2 ) is defined as 
cov{*(ti),*(t2)} = EMtOXfe)}- E{JT(t 1 )}E{JT(t 2 )} 

= <£{ (jf^j-EtJTft!)}) (x(t 2 ) - £{X(t 2 )}) }. 

(A.44) 
Parameters JT(ti) and X(t 2 ) are said to be uncorrelated if their covariance is van- 
ished, or equivalently 



var{j:(t 1 )+X(t 2 )} = var{X(t 1 )} + var{X(t 2 )}; 



(A4.5) 



otherwise, they are correlated random variables. 

When the covariance function of two random variables is normalised by their 
variances, the outcome is named the correlation coefficient and is defined as 



e{jc(ti),jc(t 2 )} 



cov{jr(ti),jr(t 2 )} 

[var{^(t 1 )}var[Jr(t2)}] 



(A4.6) 
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A. 4 . Stochastic Processes 
A.4.1.3 Auto-Covariance and Auto-Correlation 



Applying Equation (A.4.4) to a stationary process, the covariance of two random 



variables with lag t is given by 

v{X{t),X(t + x)}='E{ (x(t)-u(t)) (x(t + T)-u(t + T)] }. (A.4.7) 



COV1 



It can be shown that this covariance is a only a function of t. Therefore, it can 
be written as 

^.(T)=COv{Jt(t),Jt(t+T)}. (A.4.8) 

The function ^ (t) is known as the auto-covariance of the process {^"(t)}. 

The normalised value of %^ (t) is named the auto-correlation function and is 
defined as 

PM = ^-y (A.4-9) 

which may be also written in terms of the correlation coefficient as 

p(T} = e{Jr(t),*(t + T)}. (A.4.10) 

In the case of a stationary process with complex random variables, the associ- 
ated auto-covariance function is defined by 

£(t) = cov{X(t),X{t + r)} 

(A.4.11) 
= £{ (*(t) - £{*(t)}) * [jT(t + T)-E{JT(t + T)}j }, 

where the sign ( *) denotes the complex conjugate. 

A.4.2 Stationary Bivariate Processes 

So far, the statistical properties of stationary processes in their univariate form 
have been presented. We now focus on the relationship between two stationary 
processes. 

Consider two stationary univariate processes with continuous complex ran- 
dom variables, {JCi(t)} and {X 2 (t)}. The process {Jfi(t), Azft)} is known as a sta- 
tionary bivariate process if cov{Xi(t), Xj[t + t)} is a pure function of t. 
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In a similar way to equation Equation (A.4.11), the auto-covariance functions 
of {Xi(t)} and {X 2 (t)} can be defined as 

£h(t) = cov{JCi(t),JCi(t + T)}, (A.4.i2a) 



%. 22 {x) = cov{X 2 (t),X 2 (t + t)}; (A. 4 .i2b) 

respectively. The corresponding auto-correlation functions are given by 

A. lll u J 

p 22 (x) = f^f (A.4..,b) 

A.4.2.1 Cross-Covariance and Cross-Correlation 

To describe the correlation 'between' two distinct processes, we introduce the 
concept of cross-covariance with lag t defined by 

£ 21 (t) =cov{X 1 (t),X 2 (t + T)} 

(A.4.14) 

= £{ [jri(t)-E{jri(t)}]* [x 2 (t + T)-E{x 2 (t + T)}j }. 

The normalised value of ^ 21 M is called the cross-correlation function and is writ- 
ten 

P2i(T) = ^^ n (A.4.15) 

(^„(0)^ 22 (0)) 2 

Similar expressions exist for 3^ n an d P12 and it can be shown that 

a\i 2 M=^, 21 (-t). (A.4.16) 

A.4.3 Stationary Multivariate Processes 

All of the concepts presented for a stationary bivariate process can be generalised 
to a multivariate process. 
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A. 4 . Stochastic Processes 



A multivariate process with n sample space elements is expressed by 

{X(ti);X(t 2 );...;X(t n )} ) 

or in a more compact form by {X(t)}; where X is a column vector, containing 
the m variants of the process, defined for any t as 

/ *i(t) \ 



Xftl 



(A.4.17) 



\ x m {t) J 
Using this definition, the covariance matrix, R, at lag t may be written as 

R(T)={^ ij (T)}, i= l,2,...,m, j = l,2,...,m; (A.4.18) 

where each matrix element is given by 
£ t j(T) = cov{JC,(t),JC i (t + T)} 



£{ (x, (t) - £{Xj (t)}] * ( Jd(t + T) - £{^(t + t)}] }. 



(A.4.19) 
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Nomenclature 



Roman Letters 




& 


Complex field amplitude of light 


% 


Adiabatic bulk modulus 


c 


Thermodynamic speed of sound of fluid 


C 


Density of particle-source number 


c D 


Drag Coefficient 


c 


Speed of light 


Cp 


Specific heat of fluid at constant pressure 


c v 


Specific heat of fluid at constant volume 


7 


Displacement field 


D 


Diameter 


di 


Image distance 


d 


Displacement field in image space 


D L 


Diameter of lens 


do 


Object distance 


£ 


Event 


g 


Experiments 


!F 


Distribution function 


F 


Force 


f 


Frequency 


f# 


Lens f-number 


/ 


Probability distribution function (PDF) 


Fl 


In-plane loss-of-correlation 
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Fo Out-of-plane loss-of-correlation 

f r PDF of size distribution of aerosol particles 

f v PDF of velocity distribution of aerosol particles 

G Particle-source pattern 

G Particle-image pattern 

3i Impulse response 

3i% Impulse response of imaging system 

J£i a i Discrete impulse response 

3i G Impulse response of particle-source pattern 

I Windowed irradiance 

Axj = (x I; yj), Interrogation area dimensions in image plane 

I A Discrete irradiance field 

l A Mean of Discrete irradiance field 

3 Irradiance 

3 Intensity of light field 

j7 n The nth order Bessel function of the first kind. 

k Wave number 

k Complex wave number 

k Wave number of light field 

L Length 

At Resolution limit of lens 

I Displacement amplitude of the sound generating device 

M Magnification 

N Interrogation area size in pixels 

N F The Fresnel number 

IN" i Image density 

K p Number of particle-sources 

K p Instantaneous number of aerosol particles 

H p Instantaneous number density of aerosol particles 

K Po Equilibrium number of aerosol particles 

n po Equilibrium number density of aerosol particles 
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5^"p 


Fluctuation of number of aerosol particles 


6n p 


Fluctuation of number density of aerosol particles 


^Re 


Reynolds number of the aerosol particles 


K s 


Source density 


n s 


Scanning rate 


V 


Probability 


-Po 


Equilibrium pressure of fluid 


y 


Instantaneous pressure of fluid 


V 


Acoustic pressure 


V 


Reference pressure level 


Vo 


Acoustic pressure amplitude 


Apx 


Dimensions of pixel area 


a 
^ 


Estimated covariance matrix element 


r 


Radius 


?A 


Radius of the Airy disk 


Ar 


Resolution limit of recording medium 


1"I 


Radius of interrogation spot 


r p 


Radius of source-particle 


f p 


Radius of particle- image 


S 


Area 


s 


Correlation plane 


s 


Condensation of fluid 


^A 


= (^a,.^, )> Discrete correlation plane 


Si 


Interrogation area 


s 


Correlation plane in image space 


t 


Time 


T 


Temperature in degrees Celsius 


Tk 


Temperature in Kelvin 


tv 


Velocimetry time 


V 


Volume 


V* 


= (v X) Vy > v z ) , Velocity vector, Particle velocity of a fluid element 
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vo 


Acoustic velocity amplitude 


Vi 


Interrogation cell 


^p 


Velocity of source-particle, velocity of aerosol particles 


Vr 


Resolution cell 


w 


Window function 


Wp 


Bandwidth 


x* 


= [x,y, z), Position vector 


*,{U 


Position normal vectors 


" A 


= (x-a, Va), Discrete position vector 


X 


= (x,tj), Position vector in image space 


■t 


Random variable parameter 


X 


Random variable 


X 


Random variable vector 


Ax L 


Length of light sheet 


Z 


Acoustic impedance 


Z 


Characteristic impedance 


Az L 


Thickness of light sheet 


Z M 


Mechanical impedance 


Zs 


Specific impedance 


Greek Letters 




a 


Absorption coefficient 


<*s 


Spatial absorption coefficient 


Y 


Specific heat ratio 


r 


Valid-data yield 


Yf 


Focal length 


6 


Dirac delta-function 


11 


Dynamic viscosity coefficient of fluid 


MB 


Bulk viscosity coefficient of fluid 


d 


The angle between the light sheet normal and the tube 


K 


Thermal conductivity of fluid 


A 


Wavelength 
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A Wavelength of light field 

y Kinematic (shear) viscosity coefficient of fluid 

u Mean 

u-g Flow-around coefficient 

Up Entrainment coefficient 

Po Equilibrium density of fluid 

p Instantaneous density of fluid 

p p Density of source-particle 

p Autocorrelation 

T; Illumination duration 

t s Scanning duration 

6t s Scanning step 

a v Width of velocity distribution of aerosol particles 

T f Relaxation time of viscous fluid elements 

T p Relaxation time of aerosol particles 

t s Spatial relaxation time of viscous fluid elements 

cp Velocity potential of sound wave 

4>o Velocity potential amplitude of sound wave 

cp Phase 

cp Phase shift angle 

cp Phase of light field wave 

X = (x%. Xy )> Position vector in the Fourier's space 

X Thermal diffusivity of fluid 

£, = (£, x , E, v , E, z ), Displacement from equilibrium of a fluid element 

CI Sample space 

CI Element of sample space 

cu Angular frequency 

cu Angular frequency of light field wave 

Other Symbols 

i Imaginary unit 

• Correlation 
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© Convolution 

Scalar product 

->- Vector symbol 

a Variance 

G{ , } Correlation coefficient operator 

cov{ , } Covariance operator 

M. Covariance matrix 

% Covariance matrix element 

d Differentiation 

3 X Partial differentiation 

V Gradient 

V. Divergent 

Vx Curl 

V 2 = 5 2 + dy + 5 2 , Laplacian operator 

D 2 d'Alembertian operator 

□j-_ Lossy d'Alembertian operator 

E{ } Expectation operator 

( } Mean operator 

(} x Mean value over x 

5 [] Fourier transform operator 

2$I [] Discrete Fourier transform (DFT) operator 

Dgl -1 [] Reverse discrete Fourier transform (DFT) operator 

m [] Real part 

3[] Imaginary part 
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About the Cover 

A PIV analysis of two woodcut engravings corresponding to 1666 (left) and 1999 
(right) based on the novel The Club Dumas by Arturo Perez-Reverte. The 
1666 engraving belongs to the book De Umbrarum Regis Novum Portis written by 
a fictional author named Aristide Torchia in Venice, where the actual engrav- 
ing was performed by the artist Francisco Sole. The 1999 engraving is from the 
film The Ninth Gate by Roman Polanski. 

Aligned by the engraved borders, the images were analysed using correlation- 
based correction combined with the symmetrical phase-only filter. 
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